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Abstract 

The proton spin problem triggered by the EMC experiment and its present status 
are closely examined. Recent experimental and theoretical progresses and their im- 
plications are reviewed. It is pointed out that the sign of the sea-quark polarization 
generated perturbatively by hard gluons via the anomaly mechanism is predictable: 
It is negative if the gluon spin component is positive. We stress that the polarized 
nucleon structure function gi{x) is independent of the A;_l -factorization scheme chosen 
in defining the quark spin density and the hard photon-gluon scattering cross section. 
Consequently, the anomalous gluon and sea-quark interpretations for Fi, the first mo- 
ment of are equivalent. It is the axial anomaly that accounts for the observed 
suppression of T\. 
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1 Introduction 



Experiments on polarized deep inelastic lepton-nucleon scattering started in the middle 70s 
Measurements of cross section differences with the longitudinally polarized lepton beam 
and nucleon target determine the polarized nucleon structure function gi{x). In 1983 the first 
moment of the proton spin structure function, = Jq g^{x)dx, was obtained by the SLAC- 
Yale group to be 0.17±0.05 0, which is in nice agreement with the prediction 0.171±0.006 
expected from the Ellis- Jaffe sum rule based on the assumption of vanishing strange-sea 
polarization, i.e.. As = 0. Therefore, the polarized DIS data can be explained simply by the 
valence quark spin components. Around the same time, two theoretical analyses by Lam 
and Li [|] and by Ratcliffe were devoted to studying the gluon effects on the polarized 
proton structure function and its first moment. It appears that there is an anomalous gluon 
contribution to in the sense that, though the gluon effect is formally a leading-order QCD 
correction, it does not vanish in the asymptotic limit — ^ oo @] . However, the implication 
of this observation was not clear at that time. 

The 1987 EMC experiment then came to a surprise. The result published in 1988 and 
later indicated that = 0.126 ± 0.018, substantially lower than the expectation from the 
Ellis- Jaffe conjecture. This led to the stunning implication that very little ( < 15%) of the 
proton spin is carried by the quarks, contrary to the naive quark-model picture. While the 
proton spin arises entirely from the quarks in the non-relativistic constituent quark model, 
the sum of the z-component of quark spins, AS, accounts for 3/4 of the proton spin in the 
relativistic quark model. The EMC data implied a substantial sea-quark polarization in the 
region x < 0.1, a range not probed by earlier SLAC experiments. The question of what 
is the proton spin content triggered by the EMC experiment has stimulated a great deal of 
interest in understanding the origin of the so-called (although not pertinently) "proton spin 
crisis" . Up to date, there are over a thousand published papers connected to this and related 
topics. 

During the period of 1988-1993, theorists tried hard to resolve the proton spin enigma 
and seek explanations for the EMC measurement of T'l, assuming the validity of the EMC 
data at small x (0.01 < x < 0.1) and of the extrapolation procedure to the unmeasured small 
X region {x < 0.01). One of the main theoretical problems is that hard gluons cannot induce 
sea polarization perturbatively for massless quarks due to helicity conservation. Hence, it is 
difficult to accommodate a large strange-sea polarization As ~ —0.10 in the naive parton 
model. Right after the EMC experiment, the effect of anomalous gluon contributions to 

was revived by Efremov and Teryaev 0, Altarelli and Ross P], Carlitz, Collins and 
Mueller (see also Leader and Anselmino [|lOl)- Roughly speaking, the anomalous gluon 
effect originating from the axial anomaly mimics the role of sea polarization if the gluon spin 
component AG has a sign opposite to As. Consequently, AS is not necessarily small whereas 
As is not necessarily large. This anomalous mechanism thus provides a plausible and simple 
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solution to the proton spin puzzle: It explains the suppression of observed by EMC and 
brings the improved parton model close to what expected from the quark model, provided 
that AG is positive and large enough. But then we face a dilemma. According to the OPE 
approach, which is a model-independent approach, Ti does not receive contributions from 
hard gluons because only quark operators contribute to the first moment of at the twist-2 
and spin-1 level. This conflict between the anomalous gluon interpretation and the sea-quark 
explanation of has been under hot debate over past many years. 

In spite of much controversy over the aforementioned issue, this dispute was already re- 
solved in 1990 by Bodwin and Qiu (see also Manohar |]12[). They emphasized that the size 
of the hard-gluonic contribution to is purely a matter of the /c^-factorization convention 
chosen in defining the quark spin density Aq{x) and the hard cross section for photon-gluon 
scattering. As a result, the above-mentioned two different popular interpretations, corre- 
sponding to chiral-invariant and gauge-invariant factorization schemes respectively, are on 
the same footing. Their equivalence will be shown in details in Sec. 4 in the framework 
of perturbative QCD. The axial anomaly that breaks chiral symmetry can generate neg- 
ative helicity even for massless sea quarks. Therefore, a sizeable strange-sea polarization 
As ~ —0.10 is no longer a problem in the sea-quark interpretation. Despite this clarifica- 
tion, some of recent articles and reviews are still biased towards or against one of the two 
popular explanations for T^; this is considerably unfortunate and annoying. 

One can imagine that after a certain point it is difficult to make further theoretical or 
phenomenological progress without new experimental inputs. Fortunately, this situation 
was dramatically changed after 1993. Since then many new experiments using different 
targets have been carried out. In 1994 SMC and SLAC-E143 have reported independent 
measurements of gi{x) and confirmed the validity of the EMC data. The new world average 
AE ~ 0.30 indicates that the proton spin problem becomes less severe than before. The 
new measurements of polarized neutron and deuteron structure functions by SMC, SLAC- 
E142 and SLAC-E143 allowed one to perform a serious test on the Bjorken sum rule. This 
year marks the 30th anniversary of this well-known sum rule. We learned that it has been 
tested to an accuracy of 10% level. Data on the transverse spin structure function g2{x) 
have just become available. A probe of g2 might provide a first evidence of higher-twist 
effects. Finally, the x-dependent spin distributions for u and d valence quarks and for non- 
strange sea quarks have been determined for the first time by measuring semi-inclusive spin 
asymmetry for positively and negatively charged hadrons from polarized DIS. In short, the 
experimental progress in past few years is quite remarkable. 

On the theoretical side, there are also several fascinating progresses. For example, two 
successful first-principles calculations of the quark spin contents based on lattice QCD were 
published last year. The calculation revealed that sea-quark polarization arises from the 
disconnected insertion and is empirically SU(3)-fiavor symmetric. This implies that the axial 
anomaly is responsible for a substantial part of sea polarization. The lattice calculation also 
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suggests that the conventional practice of decomposing Aq into valence and sea components 
is not complete; the effect of "cloud" quarks should be taken into account. Other theoretical 
progresses will be surveyed in Sec. 6. 

With the accumulated data of gi{x), gi{x) and gf{x) and with the polarized two- loop 
splitting functions available very recently, it became possible to carry out a full and consistent 
next-to-leading order analysis of gi{x, Q^) data. The goal is to determine the spin-dependent 
parton distributions from DIS experiments as much as we can, especially for sea quarks and 
gluons. 

There are several topics not discussed in this review. The transverse spin structure 
function g2{x) is not touched upon except for a brief discussion on the Wandzura-Wilczek 
relation in Sec. 4.2. The small or very small x behavior of parton spin densities and polarized 
structure functions is skipped in this article. Perspective of polarized hadron colliders and ep 
colliders will not be discussed here. Some of the topics can be found in a number of excellent 
reviews [12-24] on polarized structure functions and the proton spin problem. 



2 Polarized Deep Inelastic Scattering 
2.1 Experimental progress 

Before 1993 it took averagely 5 years to carry out a new polarized DIS experiment (see 
Table I). This situation was dramatically changed after 1993. Many new experiments mea- 
suring the nucleon and dcuteron spin-dependent structure functions became available. The 
experimental progress is certainly quite remarkable in the past few years. 

In the laboratory frame the differential cross section for the polarized lepton-nucleon 
scattering has the form 

£a 1 E' , , 

-I^L,^W^\ (2.1) 



dVLdE' 2M g4 e 

where E [E') is the energy of the incoming (outgoing) lepton, L^i, and W^,^ are the leptonic 
and hadronic tensor, respectively. The most general expression of W^,^ reads 

= + M.) + - L - P-l^ I /(p . ,) (2.2) 



+ ? W'^? {s gi + 

p-q I 



p-q 




that is, it is governed by two spin- averaged structure functions Fi and F2 and two spin- 
dependent structure functions gi and 512- 
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Experimentally, the polarized structure functions gi and g2 are determined by measuring 
two asymmetries: 

da'^^ — dcr^^ da^'^ — da^'^ 

^1' " dan + d^TT ' = d^I^Td^' ^ ^ ^ 

where da^ (da^) is the differential cross section for the longitudinal lepton spin parallel 
(antiparallel) to the longitudinal nucleon spin, and da^'^ {da^'^) is the differential cross 
section for the lepton spin antiparallel (parallel) to the lepton momentum and nucleon spin 
direction transverse to the lepton momentum and towards the direction of the scattered 
lepton. It is convenient to recast the measured asymmetries A\\ and A± in terms of the 
asymmetries Ai and A2 in the virtual photon-nucleon scattering: 

A. = '-111^, A, = (2.4) 

cri/2 + cr3/2 (7-1/2 + 0-3/2 

where ai/2 and a3/2 are the virtual photon absorption cross sections for 7*(1) + iV(— ^) 
and 7*(1) + N{-^) scatterings, respectively, and a^^ is the cross section for the interference 
between transverse and longitudinal virtual photon-nucleon scatterings. The asymmetries 
Ai and A2 satisfy the bounds 

\Ai\ < 1, IA2I < VR, (2.5) 

where R = ul/ctt and (Jt = (c-i/2 + C3/2)/2. The relations between the asymmetries Ay, A± 
and Ai, A2 are given by 

A|l = D{A, + riA2), A^ = D{A2 - CA), (2.6) 

where D is a depolarization factor of the virtual photon, rj and ^ depend only on kinematic 
variables. The asymmetries Ai and A2 in the virtual photon-nucleon scattering are related 
to the polarized structure functions gi and g2 via 

A, = ^LTf*. A, = + . (2.7) 

where 'j = Q/u = Q / {E — E') = IMxj \fQ^ . Note that the more familiar relation A\ = g\jF\ 
is valid only when 7 or (?2 ~ 0. By solving (2.6) and (2.7), one obtains expressions of gx 
and g2 in terms of the measured asymmetries A\\ and A^. Since 7 — >■ in the Bjorken limit, 
it is easily seen that to a good approximation, Ay ~ DA\ and 

gx{x, Q ) _ F,{x, Q ) - - ^^^^qr^^^^g^yy ^- (2-8) 

Some experimental results on the polarized structure function gi{x) of the proton, gi{x) 
of the neutron, and gf{x) of the deuteron are summarized in Table I. The spin-dependent 
distributions for various targets are related by 

g',{x)+g-{x)= ^-^gtix), (2.9) 
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where ujd = 0.058 is the probabihty that the deuteron is in a D state. Since experimental 
measurements only cover a limited kinematic range, an extrapolation to unmeasured x 
and X — s> 1 regions is necessary. At small x, a Regge behavior gi{x) oc x"*-"-* with the 
intercept value < a(0) < 0.5 is conventionally assumed. In the EMC experiment 
a(0) = is chosen so that gi{x) approximates a constant ~ 0.2 as x < 0.01, and hence 



/o 9i{x)EMcdx = 0.002 . However, the SMC data p6| of show a tendency to rise at low 
X {x < 0.02), and it will approach a constant 1.34 ± 0.62 as x < 0.003 if a(0) = is chosen. 
Then we have Jq '^^'^ gi{x)dx = 0.004 ± 0.002 . Using the SMC data at small x and the above 
extrapolation yields Sq '^^ gi{x)sy[cdx = 0.017 ± 0.006. This explains why Ti obtained by 
SMC is larger than that of EMC (see Table I). 

Table I. Experiments on the polarized structure functions g^{x,Q'^), gi{x,Q'^) and gf{x,Q'^). 



Experiment 


Year 


Target 


m 

(GeV2) 


X range 


ptargct 

= t,gr^''\x,{Q^))dx 


Reference 


E80/E130 


1976/1983 


P 


- 5 


O.K x < 0.7 


0.17±0.05* 


m 


EMC 


1987 


P 


10.7 


0.01 < a; < 0.7 


0.126 ± 0.010 ±0.015t 


§ 


SMC 


1993 


d 


4.6 


0.006 < X < 0.6 


0.023 ±0.020 ±0.015 


n 


SMC 


1994 


P 


10 


0.003 < x < 0.7 


0.136 ±0.011 ±0.011 




SMC 


1995 


d 


10 


0.003 < a; < 0.7 


0.034 ±0.009 ±0.006 




E142 


1993 


n 


2 


0.03 < a; < 0.6 


-0.022 ±0.011 


ii 


E143 


1994 


P 


3 


0.03 < a; < 0.8 


0.127 ±0.004 ±0.010 




E143 


1995 


d 


3 


0.03 < a; < 0.8 


0.042 ±0.003 ±0.004 





* Obtained by assuming a Regge behavior Ai cx x^'^'^ for small x. 

t Combined result of E80, E130 and EMC data. The EMC data alone give T\ = 0.123 ± 0.013 ± 0.019 . 

A serious test of the Bjorken sum rule for the difference T\ — F" [Fi being defined in 
(2.13)], which is a rigorous consequence of QCD, became possible since 1993. The current 
experimental results are 



SMC [27] : F?-F^ = 0.199 ±0.038 at = iq GeV^ 
E143 H : V\- ri = 0.163 ± 0.010 ± 0.016 at = 3 GeV^ 

to be compared with the predictions!^ 

Ff-F^ = 0.187 ±0.003 with a,(10 GeV^) = 0.24 ± 0.03 , 
F?-F^ = 0.171 ±0.008 with a,(3GeV2) = 0.35 ±0.05, 



(2.10) 



obtained from the Bjorken sum rule evaluated up to for three light flavors |^ 

2 



1 flu 
6 gv 



a^m 43 ( a,m 
12 



20.22 



TT 




(2.11) 



(2.12) 



^The theoretical value 0.187 ± 0.003 for T{ - T'/ quoted by the SMC paper ^ seems to be obtained for 
three quark flavors rather than for four flavors. 
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where use of Qa/ Qv = F + D = 1.2573 ± 0.0028 has been made. Therefore, the Bjorken 
sum rule has been confirmed by data to an accuracy of 10% level. 

Recently, data on the transverse spin structure function g2 have just become available 
|]33| , |3^ . A probe of g2 might provide a first evidence of higher- twist effects. Finally, the x- 
dependent spin distributions for u and d valence quarks and for non-strange sea quarks have 
been determined for the first time by measuring semi-inclusive spin asymmetry for positively 
and negatively charged hadrons from polarized DIS [^]. For some discussions, see Sec. 8. 



2.2 The proton spin crisis 

From the parton-model analysis in Sec. 3 or from the OPE approach in Sec. 4, the first 
moment of the polarized proton structure function 

r?(Q^)= [\f{x,Q')dx, (2.13) 
Jo 

can be related to the combinations of the quark spin components via0 

r? = WelAqiQ') = E ^'(pI W5?b)^^ (2-14) 

^ q ^ q 

where Aq represents the net helicity of the quark flavor q along the direction of the proton 
spin in the infinite momentum frame: 

Aq = / Aq{x)dx = / q\x) + <f {x) — {x) — q^ {x) dx. (2.15) 

For a definition of Aq in the laboratory frame, see Sec. 4.5. At the EMC energies {Q"^) = 
10.7 GeV^ or smaller, only three light flavors are relevant 

r?(Q2) = 1 (^A«(Q^) + lAdiQ') + Ias{Q')) . (2.16) 

Other information on the quark polarization is available from the nucleon axial coupling 
constants g\ and g%: 

g\{Q^) = Au{Q^) - Ad{Q^), g\{Q^) = Au{Q^) + AdiQ"") - 2As(g2). (2.17) 

Since there is no anomalous dimension associated with the axial- vector currents A^, and Af,, 
the non-singlet couplings g\ and g\ do not evolve with and hence can be determined at 
= from low-energy neutron and hyperon beta decays. Under SU(3)-flavor symmetry, 
the non-singlet couplings are related to the SU(3) parameters F and D by 

g\=F + D, g\=3F-D. (2.18) 



^ As will be discussed at length in Sec. 4, whether or not gluons contribute to depends on the 
factorization convention chosen in defining the quark spin density Aq{x). Eq.(2.14) is valid in the gauge- 
invariant factorization scheme. However, gluons are allowed to contribute to gi{x) and to the proton spin, 
irrespective of the prescription of fcj^-factorization. 
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We use the values 



■361 



F = 0.459 ±0.008, D = 0.798 ± 0.008 , F//} = 0.575 ± 0.016 (2.19) 

to obtain = 0.579 ± 0.025 . 

Prior to the EMC measurement of polarized structure functions, a prediction for 
was made based on the assumption that the strange sea in the nucleon is unpolarized, i.e.. 
As = 0. It follows from (2.16) and (2.17) that 

r?(g^) = + ^9\- (2.20) 

This is the Ellis- Jaffe sum rule 0. It is evident that the measured results of EMC, SMC 
and E143 for Fj* (see Table I) are smaller than what expected from the Ellis- Jaffe sum 
rule: F? = 0.185 ± 0.003 without QCD corrections (= 0.171 ± 0.006 at = 10 GeV^ to 
leading-order corrections). 

To discuss QCD corrections, it is convenient to recast (2.16) to 



vrm = c^sm (^±^^?i + ^/a) + g csmgim. (2.21) 

where the isosinglet coupling is related to the quark spin sum: 

(?^(Q2) = AS(g2) = Au{Q^) + Ad{Q^) + As{Q^). (2.22) 

Perturbative QCD corrections to Fi have been calculated to 0{al) for the non-singlet coef- 
ficient Cns and to C(a^) for the singlet coefficient Cs [0, |37|:P| 

„ ^„2^ ' '^■■^ 43 /a, \^ /a,^^ 

TT iZ V TT / V TT 

\ 2 

12' 



1 - 


as 


43 


TT 


12 V TT j 


1 - 


«s 
TT 


1.10 

\ TT 


= a. 


iQ')- 





CsiQn = 1-^-1.10^ (2.23) 



From (2.17), (2.18) and the leading-order QCD correction to F^ in (2.21) together with 
the EMC resuh F?(Q2) = 0.126 ± 0.010 ± 0.015 g, we obtain 

Am = 0.77 ±0.06, Arf = -0.49 ± 0.06 , As = -0.15 ± 0.06 , (2.24) 



■^The singlet coefHcient is sometimes written as 

1 a.s 



6 TT V TT / 



in the literature, but this is referred to the quark polarization in the asymptotic limit, namely Aq{Q^) 
Aq{(X)). The above singlet coefficient is obtained by substituting the relation 

<?^(Q^)= + + 1.21 (-)' + •• -l ff^(oo) 

[ 6 TT V TT / J 

into (2.21). 
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and 



AS = 0.14 ±0.17 (2.25) 

at = 10.7 GeV^. The results (2.24) and (2.25) exhibit two surprising features: The 
strange-sea polarization is sizeable and negative, and the total contribution of quark helicities 
to the proton spin is small and consistent with zero. This is sometimes referred to as (though 
not pertinently) the "proton spin crisis" . 

The new data of SMC, E142 and E143 obtained from different targets are consistent with 
each other and with the EMC data when higher-order corrections in (2.21) are taken into 



account |3^. A global fit to all available data evaluated at a common Q in a consistent 



treatment of higher-order perturbative QCD effects yields |^ 



Am = 0.83 ±0.03, Arf = -0.43 ± 0.03 , As = -0.10 ± 0.03 , (2.26) 

and 

AE = 0.31 ±0.07 (2.27) 
at = 10 GeV^. An updated analysis with most recent data (mainly the E142 data) gives 



Am= 0.82±0.03±--- , Ad = -0.44 ± 0.03 ± ■■■ , As = -0.11 ± 0.03 ± (2.28) 
and 

AS= 0.27±0.04±--- (2.29) 

at = 3GeV^, where dots in (2.28) and (2.29) represent further theoretical and systemat- 
ical errors remained to be assigned. Evidently, the proton spin problem becomes less severe 
than before. Note that all above results for Aq and AS are extracted from data based on 
the assumption of SU(3)-flavor symmetry. It has been advocated that SU(3) breaking will 
leave AS essentially intact but reduce As substantially ||39|. However, recent lattice calcu- 



lations indicate that not only sea polarization is of order —0.10 but also it is consistent with 
SU(3)-flavor symmetry within errors (see Sec. 6.1). It is also worth remarking that elastic 
up scattering, which has been suggested to measure the strange-sea polarization, actually 
measures the scale-independent combination (As — Ac) instead of the scale-dependent As 
(see Sec. 8). 

The conclusions that only a small fraction of the proton spin is carried by the quarks 
and that the polarization of sea quarks is negative and substantial lead to some puzzles, for 
example, where does the proton get its spin ? why is that the total quark spin component 
is small ? and why is the sea polarized ? The proton spin problem emerges in the sense 
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that experimental results are in contradiction to the naive quark-model's picture. The non- 
relativistic SU(6) constituent quark model predicts that A-u = |, Ac? = — | and hence 
AE = 1, but its prediction = | is too large compared to the measured value 1.2573±0.0028 
[]32| . In the relativistic quark model the proton is no longer a low-lying S'-wave state since 



the quark orbital angular momentum is nonvanishing due to the presence of quark transverse 
momentum in the lower component of the Dirac spinor. The quark polarizations Am and 
Ad will be reduced by the same factor of | to 1 and — |, respectively, if g\ is reduced from 
I to I (see also Sec. 6.1) The reduction of the total quark spin AE from unity to 0.75 by 
relativistic effects is shifted to the orbital component of the proton spin so that the spin sum 
rule now reads 



i= ^iAu + Ad) + Ll. (2.30) 

Hence, it is expected in the relativistic constituent quark model that 3/4 of the proton spin 
arises from the quarks and the rest of the proton spin is accounted for by the quark orbital 
angular momentum. 

Let Ag be decomposed into valence and sea components, Aq = Aq^ + Aqg. The exper- 
imental fact that AE ~ 0.30, much smaller than the quark-model expectation 0.75, can be 
understood as a consequence of negative sea polarization: 

AE = AE^ + AE, = {Au^ + Ad^) + (Am, + A4 + As). (2.31) 

Nevertheless, we still encounter several problems. First, in the absence of sea polarization, 
we find from (2.17) and (2.18) that Auy = 0.92, A^ = -0.34 and AE„ = 0.58. As first 
noticed by Sehgal [0, even if sea polarization vanishes, a substantial part of the proton 
spin does not arise from the quark spin components. In fact, the Ellis- Jaffe prediction 
r^(lOGeV^) = 0.171 ± 0.006 is based on the above "canonical" values for Ag^, and AE^. 
Our question is why the canonical AE^, still deviates significantly from the relativistic quark 
model expectation 0.75 . A solution to this puzzle will be discussed in Sec. 6.1. It turns out 
that the canonical valence quark polarization is actually a combination of "cloud-quark" and 
truly valence-quark spin components. Second, in the presence of sea-quark polarization, the 
spin sum rule must be modified to include all possible contributions to the proton spin:0 

^ = i(An + Ad + As) + AG + LI + L^, (2.32) 

where AG = Jq AG{x)dx = /q^[G^(x) — G^{x)\dx is the gluon net hehcity along the proton 
spin direction, and L\^'^'^ is the quark (gluon) orbital angular momentum. It is a most great 
challenge, both experimentally and theoretically, to probe and understand each proton spin 
content. 



^It has been argued that in the double limit, ruq and Nc oo, where iriq and Nc are the light quark 
mass and the number of colors respectively, one has g'\{= AS) — and AG = Lf = 0, so that the proton 
spin is orbit in origin ^ . 
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Figure 1: Two theoretical curves for xc^{x). Tlie solid curve is obtained by fitting to the EMC Q 



and SMC |26| data at x > 0.2 with the polarized valence quark distributions given by (7.11), and 
the dotted curve arises from the valence spin distributions given by (2.33). At first sight, the latter 
curve appears to give a reasonable "eye-fit" to the data, even though its first moment is too large 
compared to the measured value. 



Before closing this section, we wish to remark that experimentally it is important to 
evaluate the first moment of gi{x) in order to ensure that the existence of sea polarization is 
inevitable. Suppose there is no spin component from sea quarks, then it is always possible 
to parametrize the valence quark spin densities, for example^ 



Am„(x) = 0.355 a;-°-^^(l-a; 
MJx) = -0.161 x-°-5^(l - 



+ I8.36x), 
x)^-6^(l + 18.36x). 



(2.33) 



in such a way that they make a reasonable "eye-fit" to the EMC and SMC |26| data of 
gi{x) even at small x (see the dotted curve in Fig. 1). One cannot tell if there is truly a 
discrepancy between theory and experiment unless is calculated and compared with data 
[(2.33) leads to 0.171 for Ti]. This example gives a nice demonstration that an eye-fit to the 
data can be quite misleading ^3|. Since the unpolarized strange-sea distribution is small 
at X > 0.2, the positivity constraint |As(x)| < s{x) implies that the data of gi{x) should 
be fully accounted for by Au^(x) and Ady{x). In Sec. 7.1 we show that a best least fit 
to gi{x) leads to a parametrization (7.11) for valence quark spin densities. The theoretical 
curve of gf without sea and gluon contributions is depicted in Fig. 1. It is clear that a 
deviation of theory from experiment for gi{x) manifests at small x where sea polarization 
starts to play an essential role. 



^This parameterization is taken from 
satisfy the first-moment constraint: Au^ 



|| except that we have made a different normalization in order to 
0.92 and = -0.34 . 
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3 Anomalous Gluon Effect in the Parton Model 

3.1 Anomalous gluon contributions from box diagrams 

We see from Section II that the polarized DIS data indicate that the fraction of the proton 
spin carried by the light quarks inside the proton is AS ^ 0.30 and the strange-quark 
polarization is As ^ —0.10 at Q"^ = lOGeV^. The question is what kind of mechanism 
can generate a sizeable and negative sea polarization. It is difficult, if not impossible, to 
accommodate a large As in the naive parton model because massless quarks and antiquarks 
created from gluons have opposite helicities owing to helicity conservation. This implies 
that sea polarization for massless quarks cannot be induced perturbatively from hard gluons, 
irrespective of gluon polarization. (Recall that our definition of Aq (2.15) includes both 
quark and antiquark contributions.) It is unlikely that the observed As comes solely from 
nonperturbative effects or from chiral-symmetry breaking due to nonvanishing quark masses. 
(We will discuss in Sec. 4.4 the possible mechanisms for producing sea polarization.) It 
was advocated by Efremov and Teryaev 0, Altarelli and Ross Q, Carlitz, Collins and 
Mueller (see also Leader and Anselmino [1^]) that the difficulty with the unexpected 
large sea polarization can be overcome by the anomalous gluon effect stemming from the 
axial anomaly, which we shall elaborate in this section. 

As an attempt to understand the polarized DIS data, we consider QCD corrections to 
the polarized proton structure function gi{x). To the next-to-leading order (NLO) of a^, the 
expression for gi{x) i^ 



where % is the number of active quark flavors, Ag(x) = q\x) + q^ {x) — q^{x) — q^{x), 
AG{x) = {x) — G^{x), and ® denotes the convolution 



There are two different types of QCD corrections in (3.1): the Afg term arising from vertex 

and self-energy corrections (corrections due to real gluon emission account for the InQ^ 

dependence of quark spin densities) and the other from polarized photon-gluon scattering 

[the last term in (3.1)]. As we shall see later, the QCD effect due to photon-gluon scattering 

is very special: Unlike the usual QCD corrections, it does not vanish in the asymptotic 

limit oo. The Afg{x) term in (3.1) depends on the regularization scheme chosen. 

Since the majority of unpolarized parton distributions is parametrized and ffited to data in 

^Beyond NLO, it is necessary to decompose the quark spin density into singlet and non-singlet compo- 
nents; see Eq.(7.6) for a most general expression of gi{x, Q^). 





(3.1) 




(3.2) 
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the MS scheme, it is natural to adopt the same regularization scheme for polarized parton 
distributions in which [QQ 



/ln(l 



X] 



\ 1 — X 



2(1 -x) 



1+x^ 

1 — X 



Inx 



+3 + 2a;- + 6(1 -x) 



1(1 +x) 



where the "+" distribution is given by 



9{x) 

V 1 - 



dx 



' fix) '-^^l^^dx. 

1 — X 



(3.3) 



(3.4) 



The first moment of fq{x) and Afq{x) is and —2, respectively. Note that the first moment 
of A/g(x) is scheme independent at least to NLO. 




k + permutations 



Figure 2: The photon-gluon scattering box graph. 

To compute the polarized photon-gluon cross section Aa'^'^ amounts to evaluating the 
box diagram of photon-gluon scattering with a physical cutoff on the intermediate states 
(see Fig. 2). Using the relation 



(3.5) 



^The expression of (3.3) for Afq{x) is valid for /il^^.^ = Q^, where fifact is a factorization scale to be 
introduced below. When /^f^^.^. ^ Q^, the contribution Q 



In 



should be added to IS.fq{x,Q'^). 
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where e^{±) = (0, 0, 1, ±i)/-\/2 is the transverse polarization of external gluons and eij is an 
antisymmetric tensor with ei2 = 1, the contribution of Fig. 2 to Acr^*^ for a single quark 
flavor is HI 



a 



qq) - a(7(+)G(-) ^ qq) 



cPk^dk+dk- 2^5[{p - kf - m^]5[{q + kf - m^]Ti{- ■ ■} 



(/c^ — rn? + iey 



(3.6) 



with 



Tr{---} = ^ei,Q„Tr{7^'(^ + m)7"(^ + |^ + m)7'(^ + m)7^(jf^-|^-m) 

+7'(l^ + ^-^ + m)-f''{^ + m)Y(^ + m)^\l^ + ^ + m)}, (3.7) 

where T = |, m is the mass of the quark, p is the gluon momentum, and is the quark 
transverse momentum perpendicular to the virtual photon direction. It is convenient to 
evaluate the integral of (3.6) in the light-front coordinate = {p~ ,p'^ ,p±) with p± = and 
p± = (pO ±p3)/y/2. a tedious but straightforward calculation yields (for a derivation, see 
e.g., for the general case and [H for on-shell gluons): 



Ztt Jo 



A-2 



'l-2x){kl+m'^) -2m^{l-x) 1 - 2x' 



kl/K^ 



[ k^ + — p'^x{l — x)]'^ 



2K^ 



where = [(1 — x)/Ax]Q'^ and x = k~^/p~^. Note that all higher-twist corrections of order 
p^ IQ"^ and rn? IQ^ have been suppressed in (3.8). It is evident that I\o'^'^[x) has infrared and 



collinear singularities at vn? 



p" 



and A;^ = 0. Hence, we have to introduce a soft cutoff 



to make Aa"''~'{x) finite. For » m^, —p"^, Eq.(3.8) reduces to (for an exact expression of 
Aa'^'^(x) after k± integration, see |^ W^ ) 



Aa-^{x,Q') = ^ 
Ztt 



(2x 



In 



x] 



rn? — p^x(l — x] 
2rn? — p^x{l — 2x) 



+ lnl^-l 

X 



m? — p^x(l — x) 
Depending on the infrared regulators, we have 



(3.9) 



AcrX^(x,Q2) 
Aaf(x,Q2) 



In + In ^ 



as 
2^ 



{2x - 1) 

(2x - 1; 



In 



-p^ 



x" 



+ In- 



X 



- 1 



In — h In 



1 — a; 



2(1 -x) 
f 2(1 - X 



(3.10) 
(3.11) 

(3.12) 



for the momentum regulator (p^ 7^ 0) the mass regulator (m^ 7^ 0) and the dimensional 
regulator [fi ^ 0) in the modified miminal subtraction scheme P], respectively. Note that 
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the coefficient {2x — 1) in Eqs.(3.10-3.12) is nothing but the spin sphtting function 2APqG{x) 
[see (3.26)] and that the term proportional to 2(1 — x) in (3.11) and (3.12) is an effect of 
chiral symmetry breaking:^ It arises from the region where k'j_ ~ in the mass-regulator 
scheme, and from k"^, ~ in the e = n — 4 dimensions in the dimensional regularization 

MS 

scheme due to the violation of the identify {7^, 75} = 0. For the first moment of Aa'^'^{x), 
it is easily seen that 

^(^ccMix)dx = ^ Aal^{x)dx = Aa][^{x)dx = 0. (3.13) 

The result (3.13) can be understood as follows. The cutoff-dependent logarithmic term, 
which is antisymmetric under x — 1 — x, makes no contribution to Jq Aa^^(x)dx, a conse- 
quence of chiral symmetry or helicity conservation. As a result, Jq Aa^'-'(x)dx receives "hard" 
contributions from ~ in the momentum-regulator scheme, but it is compensated by 
the soft part arising from k'j_ ~ in the mass-regulator scheme. 

It is clear that the cross sections given by (3.10-3.12) are not perturbative QCD reliable 
since they are sensitive to the choice of the regulator. First, there are terms depending loga- 
rithmically on the soft cutoff. Second, the first moment of Aa'^'^{x) is regulator dependent. 
It is thus important to have a reliable perturbative QCD calculation for Aa'^'~'{x) since we 
are interested in QCD corrections to gi{x). To do this, we need to introduce a factorization 
scale /ifact, so that 

Aa^^ix, Q') = Aa,t,(x, Q^ /iLt) + ^^lli^, /^Lt) (3-14) 
and the polarized photon-proton cross section is decomposed into 

i 

+ A<^,(x,Q^/iLt)®AG'(x,/iLt))• (3.15) 

That is, the hard piece of Aa'^'~'{x) which can be evaluated reliably by perturbative QCD 
contributes to gi{x), while the soft part is factorized into the nonperturbative quark spin den- 
sities Aqi{x). Since Aa'^'P{x) is a physical quantity, a different factorization scheme amounts 
to a different way of shifting the contributions between Acrh^j(x) and Aq{x). An obvious 
partition of Aa'^'^{x) is that the region where k'j_ > nf^^^ contributes to the hard cross section, 
whereas the soft part receives contributions from k\ < n}^^^ and hence can be interpreted as 
the quark and antiquark spin densities in a gluon, i.e., ActJo'^^(x, /ifa^t ^?/^fact)- Phys- 

ically, the quark and antiquark jets produced in deep inelastic scattering with k^ < nf^^^ are 
not hard enough to satisfy the jet criterion and thus should be considered as a part of one-jet 

^The 2(1 — x) term in (3.11) and (3.12) was neglected in the original work of Altarelli and Ross |^ and of 
Ratcliffe [|. One may argue that since this contribution is soft, it will not contribute to "hard" Aa'''-^ . As 
shown below, the cross sections A(t2^(x) and Aa^'^(x) without the 2(1 — x) term indeed give the correct 
result for the first moment of ActJJ'^.j in the chiral-invariant factorization scheme. 
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cross section The choice of the "ultraviolet" cutoff for soft contributions specifies the k± 
factorization convention. There are two extremes of interest: the chiral- invariant scheme in 
which the ultraviolet regulator respects chiral symmetry, and the gauge-invariant scheme in 
which gauge symmetry is respected but chiral symmetry is broken by the cutoff. 

The next task is to compute Aa]^^{x). It can be calculated from the box diagram by 
making a direct cutoff ~ /ifact on the k± integration. Note that for k'j_ < /ifa^t) box 
diagram for photon-gluon scattering is reduced under the collinear approximation for the 
quark-antiquark pair created by the gluon to a triangle diagram with the light-front cut 
vertex 7^7^ combined with a trivial photon-quark scattering As a result, Aq^{x) 

can be also obtained by calculating the triangle diagram with the constraint k'j_ < yUfact- In 
either way, one obtains 

{kl + m^){l - 2x) - 2m^{l - x) 



2n 



dki 



[k\+ m? 



-,(3.16) 



p^x{l — x)Y 

where Oil/Q"^) corrections are negligible for << and the subscript CI indicates that 
we are working in the chiral-invariant factorization scheme. The result is 



(2x - 1) In 



/^fact + ^ -P a:(l - 
m? — p'^x{l — x) 



X] 



2m? — p'^x{l — 2x) 
im? — p'^x{l — x) 



A'-fact 



X) 



(3.17) 



For /if. 



act 



Aggi(x,/i^, 



p 



'2, it reduces to 



act 7 



2-K 



2n 
2-K 



{2x-l) (ln(/iLt/-p') + ln^ 
{2x-l)H^^lJm^) + 2{l-x)]■ 
{2x-l)H^^lJ^^U + 2{l 



X] 



l-2x 



(3.18) 



for various soft cutoffs. Note that, as stressed in the soft cross sections or quark spin 
densities in a helicity + gluon given by (3.17) or (3.18) do not make sense in QCD as they 
are derived using perturbation theory in a region where it does not apply. Nevertheless, it 
is instructive to see that 



AgS'j = / Aq^j{x)dx = 0, for = 
Jo 



or — p^ >> m^. 



(3.19) 



as expected. Hence, a sea polarization for massless quarks, if any, must be produced nonper- 
turbatively or via the anomaly (see Sec. 4.4). Now it does make sense in QCD to subtract 
A^soft froni Aa'^''-^ to obtain a reliable perturbative QCD result for AcxJJ^j: 



A<frd(a;,Q^/^Lt)cI = ^(2x 

Ztt 



A'-fact 



X 



(3.20) 



Evidently, A(Th^d(a;) is independent of the infrared regulators as long as fif^^^^ » m^, —p"^; 
terms depending on soft cutoffs are absorbed into the quark spin densities. It is also clear 



17 



that the soft 2(1 — x) term in (3.11) and (3.12) drops out in Acr^^^(x). Therefore, 



hardV^'Q ) A''fact)ci — r, • 

ZTT 



(3.21) 



Since gauge invariance and hehcity conservation in the quark-gluon vertex are not broken 
in the chiral-invariant factorization scheme, it is evident that Ag^j does not evolve, con- 
sistent with the naive intuition based on hehcity conservation that the quark spin Ag^j = 
/q^ IS.cf^^{x)dx for massless quarks is independent. 
Substituting (3.3) and (3.20) into (3.1) leads to 



r?(g^) ^ dxg{{x, Q') = ^ (i - f ) E 

where Agci(Q^) = ^QciiQo) and use has been made of 



Ag.(g^)ci-^#^AG(g^) 



27r 



(3.22) 



Jo Jx y \y J Jo 



dxf{x) / dyg{y). 



(3.23) 



The (1 — ^) term in (3.22) comes from the QCD loop correction,! while the a^AG term arises 
from the box diagram of photon-gluon scattering. If the gluon polarization inside the proton 
is positive, a partial cancellation between Ag and ^AG will explain why the observed T\ is 
smaller than what naively expected from the Ellis- Jaffe sum rule. It is tempting to argue that 
the box-diagram QCD correction is negligible at large since vanishes in the asymptotic 
limit Q"^ oo. However, it is not the case. To see this, consider the Altarelli-Parisi (AP) 
equation for flavor-singlet polarized parton distribution functions: 

d_ /ASci(x,t)\ ^ «^ / ApS(x) 2nfAP,G{x)\ /ASci(x,t)\ 
dt\ AG{x,t) J 2n [aPg,{x) APgg{x) ) [ AG{x,t) )' 

where t = In^Q"^ / Aqq^j) . Although the leading-order splitting functions in 



(3.24) 



AP{x) = AP(°)(a;) + — AP(i)(x) + 

2n 



(3.25) 



have been obtained long time ago |5Q], the NLO results are not available until very recently 



||51|| . To the leading order, the AP splitting functions read 



4 / 1 + x2 
3 



1 



X 



AP^a(^^ 



V X 



Ml 

qG 
1 



\(2x 



X 



+ ^6il-x) 



3(2 -x). 



(3.26) 



with pQ = (33 - 2n/)/3. Since a^iQ^) = 47r/(/3o In QVAqcd), it follows from (3.24) that 



dlnQ^ 



{as{Q')AG{Q')) = ^AS = 0(«^ 



(3.27) 



^Recall that perturbative QCD corrections to have been calculated up to the order of a^; see Eq.(2.21). 
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Consequently, ctsAG is conserved to the leading-order QCD evolution;[^ that is, AG grows 
with InQ^, whereas is inversely proportional to InQ^. Explicitly, a solution to (3.24) 
reads 

AOm = -^AE(QS + !^|/|^(aG(QS) + ^AE(Qs). (3.28) 

Hence, hard gluons contribute to the first moment of gi{x) even in the asymptotic limit. As 
we shall see below, it is the axial anomaly that makes this QCD effect so special. 

Physically, the growth of the gluon spin with can be visualized in two different ways. 
From (3.26) we have Jq APQ^{x)dx = 2 . This means that a polarized quark is preferred to 
radiate a gluon with helicity parallel to the quark spin. Since the net quark spin component 
within the proton is positive, it is clear that AG > at least for the gluons perturbatively 
emitted from quarks. As increases, the number of gluons with + helicity radiated from 
polarized quarks also increases and this explains why AG grows with Q^. Alternatively, 
this growth also can be understood by considering the splitting of a helicity + gluon into a 
quark- antiquark pair or into two gluons. Since 

J^' AP${x)dx = 0, J^' APj^lWx = ^/3o, (3.29) 
the gluon helicity has a net gain with probability 11/2 — nj/3 > in the splitting Thus 



the gluon spin component increases with increasingly smaller distance scale. Now we see that 
perturbative QCD provides all the necessary ingredients for understanding the smallness of 
AS. As a result of anomalous gluonic contributions to Ti in the chiral- invariant factorization 
scheme, what measured in polarized DIS experiments is not Ag, but rather a combination 
of Aq and a, AG [cf. Eq.(3.22)]: 

Aq^ Aq- -^AG. (3.30) 
27r 

Consequently, (2.26) and (2.27) are replaced by 



and 



Auci ^AG = 0.83 ±0.03, 

27r 

Add ^AG = -0.43 ±0.03, (3.31) 

27r 

Asci-^AG = -0.10 ±0.03, 



(Am + Ad + As)ci - i^AG = 0.31 ± 0.07 (3.32) 



^°This constant behavior for a^AG also can be seen from the analysis of anomalous dimensions of the 
Chern-Simons current (see Sec. 4.6). 
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at g2 = lOGeV^ (3.31) and (3.32) imply that in the presence of anomalous gluon con- 
tributions, ASci is not necessarily small and Asci is not necessarily large. In the ab- 
sence of sea polarization and in the framework of perturbative QCD, it is easily seen that 
AG ~ 0.10(27r/«,) ~ 2.5 at = 10 GeV^ and AEci ~ 0.60. It thus provides a nice and 
simple solution to the proton spin problem: This improved parton picture is reconciled, to 
a large degree, with the constituent quark model and yet explains the suppression of Ti, 
provided that AG is positive and large enough. This anomalous gluon interpretation of 
the observed Ff, as first proposed in 0, |§, 0|, looks appealing and has become a popular 
explanation since 1988. Note that AG ~ 2.5 ought to be regarded as an upper limit for the 
magnitude of the gluon spin component within a proton, as it is derived by assuming no 
intrinsic strange-sea polarization (see also Sec. 4.4). 

3.2 Role of the axial anomaly 

In order to understand the origin of the anomalous gluon contribution to Ff, we consider an 
important consequence of the OPE which requires that 

f^Aa^^{x)dx= ^F+, (3.33) 

where T'^ is the contribution of the triangle diagram for the axial-vector current be- 
tween external gluons (see Fig. 3 in Sec. 4.2) evaluated in the light-front coordinate = 
{p~ ,p~^ ,p±)- The relation (3.33) ensures that the two different approaches, the OPE and 
the improved parton model, yield the same results. It has been shown in |^ that the in- 
tegrands of both sides of (3.33) are equal in the low k\ region. This in turn implies that 
AcrJo^t(x) = Ag'-^(x), namely the soft part of the photon-gluon scattering cross section equals 
to the soft part of the triangle diagram, a relation which we have employed before for com- 
puting the quark spin densities inside a polarized gluon [see Eq.(3.16)]. Moreover, we have 
shown that Aa]^i^^{x)dx = Jq Aq§^{x)dx = for = or -p^ >> [cf. Eq.(3.19)]. 
This means that only the integrands at large k\ region contribute to (3.33). 

It is well known that the triangle diagram has an axial anomaly manifested at k\ ^ oo 
(see Sec. 4.3). Since only the k'j_ ~ region contributes to the nonvanishing first moment 
of Aa'^'^{x) it follows from (3.33) that the anomalous gluon contribution — ^AG to Ff 
is intimately related to the axial anomaly. (Both sides of (3.33) have values — as/27r.) That 
is, the gluonic anomaly occurs in the box diagram (Fig. 2) at k\ = = [(1 — x)/Ax]Q'^ 
with X ^ and contributes to the first moment of Acr^^^(x). This means that the upper 
quark line in the box diagram has shrunk to a point and this point-like behavior measures 
the gluonic component of the quark Fock space Q], which is identified with the contribution 
-^AG to F?. 

At this point, it is instructive to compare unpolarized and polarized structure functions. 
The unpolarized structure function Fi(x,Q^) has a similar expression as Eq.(3.1) for the 
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polarized one. However, the first moments of unpolarized fq{x) [cf. Eq.(3.3)] and 
vanisli so tliat QCD corrections to / Fi{x) reside entirely in the evolution of the first 
moment of the unpolarized quark distributions: 



It is mainly the anomalous gluon contribution that makes / gi{x)dx behave so differently 
from / Fi{x)dx. We conclude that it is the gluonic anomaly that accounts for the disparity 
between the first moments of gi{x) and Fi{x). 

We should remind the reader that thus far in this section we have only considered the 
chiral-invariant factorization scheme in which a brute-force ultraviolet cutoff on the kj_ in- 
tegration is introduced to the soft part of the box diagram. In this case, the axial anomaly 
manifests in the hard cross section for photon-gluon scattering. However, this is not the 
only fc_|_-factorization scheme we can have. In the next section, we shall see that it is equally 
acceptable to choose a (gauge- invariant) factorization prescription in which the anomaly is 
shifted from Aa^^^^ to the quark spin density inside a gluon. Contrary to the aforemen- 
tioned anomalous gluon effects, hard gluons in the gauge-invariant scheme do not make 
contributions to the first moment of A(T2^^j(a;). 

Before ending this section we would like to make two remarks. The first one is a historical 
remark. 

1) . The first consideration of the hard gluonic contribution to was put forward by 
Lam and Li long before the EMC experiment. The questions of the regulator dependence 
in the evaluation of the photon-gluon scattering box diagram, the identification of AG with 
the forward nucleon matrix element of the Chern-Simons current (see Sec. 4.5), and the 
behavior of as AG, ■ ■ ■ etc. were already addressed by them. A calculation of Aa'^'-^{x) using 
the dimensional regularization was first made by Ratcliffe also before the EMC paper. 

2) . We see from (3.26) that Jq APqq{x) = 0. This indicates that AEci is independent. 
Physically, this is because the quark helicity is conserved by the vector coupling of a gluon to 
a massless quark. However, ASci and Agci cannot be written as a nucleon matrix element 
of a local and gauge-invariant operator (this will be discussed in Sec. 4.5). Since ASci 
does not evolve and since AG induced by gluon emissions from quarks increases with Q^, 
conservation of angular momentum requires that the growth of the gluon polarization with 

be compensated by the orbital angular momentum of the quark-gluon pair so that the 
spin sum rule (2.32) is independent; that is, Ll + L'^ also increases with with opposite 
sign. It is conjectured in Sec. 6.3 that LI in the chiral-invariant factorization scheme could 
be negative. 




(3.34) 



4 Sea Polarization Effect in the OPE Approach 
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4.1 Preamble 



We see from the last section that the anomalous gluon contribution to furnishes a simple 
and plausible solution to the proton spin problem. A positive and large gluon spin component 
will help explain the observed suppression of Ti relative to the Ellis- Jaffe conjecture and in 
the meantime leave the constituent quark model as intact as possible, e.g., AS ~ 0.60 and 
As ~ 0. However, this is by no means the end of the proton-spin story. According to the 
OPE analysis, only quark operators contribute to the first moment of Qi at the twist-2, spin-1 
level. As a consequence, hard gluons do not make contributions to Fj' in the OPE approach. 
This is in sharp conflict with the improved parton model discussed before. Presumably, the 
OPE is more trustworthy as it is a model-independent approach. So we face a dilemma here: 
On the one hand, the anomalous gluon interpretation sounds more attractive and is able 
to reconcile to a large degree with the conventional quark model; on the other hand, the 
sea-quark interpretation of Fj* relies on a solid theory of the OPE. In fact, these two popular 
explanations for the data have been under hot debate over the past years. 

Although the OPE is a first-principles approach, the sea-quark interpretation is never- 
theless subject to two serious criticisms. First, how do we accommodate a large and negative 
strange-quark polarization As ~ —0.10 within a proton ? Recall that, as we have repeatedly 
emphasized, no sea polarization for massless quarks is expected to be generated perturba- 
tively from hard gluons owing to helicity conservation. Second, the total quark spin AS in 
the OPE has an anomalous dimension first appearing at the two-loop level. This means that 
AS evolves with Q^, in contrast to the naive intuition that the quark helicity is not affected 
by the gluon emission. In the last 7 years, there are over a thousand papers triggered by the 
unexpected EMC observation. Because of the above-mentioned criticisms and because of the 
deviation of the sea-polarization explanation from the quark-model expectation, the anoma- 
lous gluon interpretation seems to be more favored in the past by many of the practitioners 
in the field. 

In this section we will point out that within the approach of the OPE it is precisely 
the axial anomaly that provides the necessary mechanism for producing a negative sea- 
quark polarization from gluons. Hence, the sea-quark interpretation of F^* is as good as the 
anomalous gluon one. In fact, we will show in the next section that these two different 
popular explanations are on the same footing; physics is independent of how we define the 
photon-gluon cross section and the quark spin density. 



4.2 A mini review of the OPE 

The approach of the operator product expansion provides relations between the moments 
of structure functions and forward matrix elements of local gauge-invariant operators (for a 
nice review, see [Q). For inclusive deep inelastic scattering, the hadronic tensor Wf^i, has 
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the expression 

W,. - d'xe'^-^ip,s\ [J,{x), Mmp,s) (4.1) 

for a nucleon state with momentum p^^ and spin s^^. Since W^i, in the DIS hmit is dominated 
by the hght-cone region where ~ (but not necessarily ~ 0), the structure of the 
current product is probed near the hght cone. In order to evaluate W/^,^, it is convenient to 
consider the time-ordered product of two currents: 

V = ^/ d*xe"^-^T{J^{x)MO)) (4.2) 

and the forward Compton scattering amplitude = (p, s), which is related to the 

hadronic tensor by the relation W^,^ = ;^ImT^,^ via the optical theorem. 

In the limit g — > oo, the operator product expansion allows us to expand t^i, in terms of 
local operators; schematically, 

1™ V = EC/'-,n(5)C'-(0)- (4-3) 

q—*oo ^ 

The Wilson coefficient functions C„ can be obtained by computing the quark or gluon matrix 
elements of J/^Ji, and 0„. Consider t^i, in the complex cu {= 1/x — 2p ■ q/ — (f) plane. By 
analyticity, the Feynman amplitude M^,^ corresponding to the free quark (or gluon) matrix 
element of J^Ji, can be calculated at uj near (but not in the physical region 1 < a; < oo) 
and expanded around uj — Q. Generically, 

M^, = (A;, A|V|A;, A) - A|0„|A;, A) (4.4) 

n 

for a quark state with momentum and spin A^. Since the free quark matrix elements of 
the quark operators have the form 

(A;,A|0(;^-''"|A;,A) = A;'^i---A;'^", (A;, A|0^^-^"|A;, A) = /lA;"! • • • A;^", (4.5) 

for vector and axial-vector types of quark operators, where /i is a helicity of the quark state, 
the coefficient functions C^j,,„ arc thus determined. 

The leading-twist (= dimension— spin) contributions to the antisymmetric (spin-dependent) 
part of t^y in terms of the operator product expansion are 

V]= E ^W-^^fA) ?Mi---?Mn-iE2C.n05''--'^"->, (4.6) 

n=l,3,--- \ y / i 

where the sum is over the leading-twist quark and gluon operators. The twist-2 quark 
and gluon operators are given by 

^.{.....,„_,} ^ 1 (4.7) 

Qa{^.i-nn-i} ^ Tr (e'^^^G^^D^^i-D^-^G^-^^) , (4.8) 
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with G^v a gluon field strength tensor, a covariant derivative, and {■ • ■} a complete 
symmetrization of the enclosed indices. The corresponding Wilson coefficients in (4.6) to 
the zeroth order of a, are 



Ci^n = e^, for quark operators of flavor g, 
C2,n = 0, for gluon operators. 



(4.9) 



It is useful to decompose the operator Oa into a totally symmetric one and a one with mixed 
symmetry 



^a{^J.l■■■^ln-l} _ /-){o-/tl--Mn-l} I /o[o"{A*l]M2---Mn-l} 

Ua — ^'^A 5 



where [■ ■ ■] indicates antisymmetrization. 



O 



{o"AilM2---Mn-i} 



n 



(4.10) 



(4.11) 



for n = 1, 3, 5, ■ ■ ■, is a twist-2 operator, and 



O 



[cr{/xi]/J2---Mn-l} 
A 



1 



n 



O 



(T{/il---/l„_l} Q^ll{a^J.2---lJ.n-l} 



(4.12) 



for n = 3, 5, ■ ■ ■, is a twist-3 operator. The proton matrix elements of these two operators 



are 



{P.s\0\[ 



WlJ.lfJ.2---fJ.n-l} 



\P, s) 



n 



(4.13) 



where aj,„ and di^n are unknown reduced matrix elements. 
Writing 

J [m^^] = «fi'l h «fi'2 7 7^ 

p ■ q [p ■ q) 

in analog to W[pi,-j [see Eq.(2.2)] and comparing with the proton matrix element of [cf. 
Eq.(4.6)] gives 



(4.14) 



91 — X/ X/ 2Cj_„aj^„ct; 

71=1, 3, 5--- i 

^ r/1 - 

^2 = 2. L 

n=l,3,5--- i 



n 



n 



n — 1 



n 



(4.15) 



It follows from (4.15) that 



du'. 



92{uj) = -gi{uj) + —gi{uj')+ 

Jo id 



n=l,3.5, 



n — 1 



n 



(4.16) 
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Using dispersion relations to relate gi^2 in the unphysical region {u ~ 0) to their values in 
the physical region {1 < uj < oo) finally yields the moment sum rules: 

/■I 1 

/ dxx"'~^gi{x) = -VCi,„aj,„, n = l,3,5, (4.17) 
Jo 2 . 

/ dx x"^"^ g2{x) = --(- ]y2Ci,n{ai.n~ di^n) , n = 3,5,---, (4.18) 

Jo 2 \ n / ^ 

and the relation]^ 

g2ix)=g^'^ix)+U^), (4.19) 

obtained from (4.16), where 

g^^(x) = -g^ix)+ (4.20) 
Jx y 



is a contribution to g2{x) fixed by gi{x), first derived by Wandzura and Wilczek ||55|, and 
g2{x) is a truly twist-3 contribution related to the twist-3 matrix elements di^n- 

For n = 1, the moment sum rule (4.17) for gi is particularly simple: Gluons do not 
contribute to the first moment of gi as it is clear from (4.8) that there is no twist-2 gauge- 
invariant local gluonic operator for n = 1, as stressed in Since 



{p,s\Ol^\p,s)= a,,,s^ (4.21) 

from (4.13) and Ci^i = to the zeroth order of [see (4.9)], it follows that 

/•II X 
gi{x)dx = -Ci,iai,i = -- (p, s| ^ 6^57^75^!^, s)s''. (4.22) 

Denoting 

{P,s\q-f^-f5q\p,s) = Sf,Aq, (4.23) 
(4.22) leads to the well-known naive parton-model result [cf. Eq.(2.16)] 

gP,{x)dx = ^ i^Au + ^Ad + ^As), (4.24) 

which is rederived here from the OPE approach. 

^^It should be stressed that the relation (4.19) is derived from (4.16) rather than from (4.17) and (4.18). 
It has been strongly claimed in that (4.19) is a priori not reliable since its derivation is based on the 
dangerous assumption that (4.17) and (4.18) are valid for all integer n. Obviously, this criticism is not 
applied to our case and (4.19) is valid as it stands. 
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4.3 Axial anomaly and sea-quark polarization 

Contrary to the improved parton model discussed in Sec. Ill, we see that there is no any 
gluonic operator contributing to the first moment of gi{x) according to the OPE analysis. 
The questions are then what is the deep reason for the absence of gluonic contributions to 
rf and how are we going to understand a large and negative strange-quark polarization ? 
The solution to these questions relies on the key observation that the hard cross section 
^^hard(^) hcucc the quark spin density Aq{x) are A;_L-factorization scheme dependent. 
We have freedom to redefine Aa^^^^{x) and Aq{x) in accord with (3.15) but the physical 
cross section Aa'^P{x) remains intact. Therefore, there must exist a factorization scheme 
that respects the OPE: Hard gluons make no contribution to and Aq can be expressed 
as a nucleon matrix element of a local gauge-invariant operator. In this scheme, gluons can 
induce a sea polarization even for massless quarks. This can be implemented as follows. As 
discussed in the last section, the quark spin density inside a gluon Aq'^{x) can be obtained 
by calculating the triangle diagram with an ultraviolet cutoff to ensure that ^ /^fact- 
is well known that in the presence of the axial anomaly in the triangle diagram, gauge 
invariance and chiral symmetry cannot coexist. So if the ultraviolet regulator respects gauge 
symmetry and axial anomaly, chiral symmetry will be broken. As a consequence, quark- 
antiquark pairs created from the gluon via the gluonic anomaly can have the same helicities 
and give rise to a nonvanishing Aq'^{x). Since the axial anomaly resides at A;^ — > oo, 
evidently we have to integrate over from to oo to achieve the axial anomaly and hence 
chiral-symmetry breaking, and then identify the ultraviolet cutoff with //fact- We see that the 
desired ultraviolet regulator must be gauge-invariant but chiral- variant owing to the presence 
of the QCD anomaly in the triangle diagram. Obviously, the dimensional and Pauli-Villars 
regularizations, which respect the axial anomaly, are suitable for our purposes. 



7 75 



P 




P 



Figure 3: The triangle graph for between external gluons. 




(4.25) 
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where T = ^, = (0, 0, 1, ±2)/-\/2 is the transverse polarization of external gluons, the 
factor 2 comes from the fact the gluon in Fig. 3 can circulate in opposite direction, and the 
dimensional regularization is employed to regulate the ultraviolet divergence. The quark 
spin density inside a gluon is then given by 

2ig^T r d'^-'^k\dk+dk- 



(2vr) 



-5 X 



k+\ Tr{---} 

J9+ / {k"^ — im? + i^Yiiv — k)'^ — m'^ + ie] 



(4.26) 



Note that Jq AqQi(x)dx = T'^ / (2p^) [cf. Eq.(3.33)]. We first perform the k integral in 
(4.26) by noting that a pole of k~ locating at 



k- 



(4.27) 



2{p+ - k+) 

in the region < A;^ < contributes to contour integration. The result for "+" helicity 
external gluons is 



x] 



as 
2^ 



X 



{kl + m^){l- 



2x) 



2m^{\ 



X] 



n 



4 



n 



kl{l-x) 



(4.28) 



where the subscript GI designates a gauge-invariant factorization scheme. The last term 
proportional to {n — 4) arises from the use of 75 in dimensional regularization. The 75 
matrix (= i7°7^7^7^) anticommutes with the Dirac matrix in 4 dimensions but commutes 
with the Dirac matrix in n — 4 dimensions. This term originating from the axial anomaly 
thus survives at k\ 00. By comparing (4.28) with (3.16), it is clear that AgQj(x) has the 
same expression as that of Ag^j(x) except for the presence of an axial-anomaly term in the 
former. It follows that 



Aqg{x) - Aqg{x) 



'27r 



{2x 



Din 



m 



p'^x{l 



XI 



A'-fact 



2 q-x) 



/^fact + "^^-P^2;(l-x)^ 



(4.29) 



for mass and momentum cutoffs, and 



XI 



XI 



'2^ 



(2x 



/^fact 



2/i 



fact(l 



X) 



A^fact + Z^— 



(4.30) 



for the dimensional infrared cutoff. Hence, 



Aq, 



X) 



Aggi(x 



1-x) 



(4.31) 



for fif^^^ » /i^, m^, — The difference between the quark spin densities in gauge- invariant 
and chiral-invariant factorization schemes thus lies in the gluonic anomaly arising at the 
region /c^ ~ /^fact- noted in passing, the quark spin distribution in a gluon cannot 
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be reliably calculated by perturbative QCD; however, the difference between AgQj(x) and 
Aq^i{x) is trustworthy in QCD. It is interesting to see from Eqs. (4.31) and (3.19) that 

Aggi(/iLt) = -^4^ (4.32) 

for massless quarks. Therefore, the sea- quark polarization perturbatively generated by helicity 
+ hard gluons via the anomaly mechanism is negative ! In other words, a polarized gluon is 
preferred to split into a quark-antiquark pair with helicities antiparallel to the gluon spin. As 
explained before, chiral-symmetry breaking induced by the gluonic anomaly is responsible 
for the sea polarization produced perturbatively by hard gluons. 

Since /^.a^^^^^x) = Aa'^'^{x) — Ag'^(x), it follows that the hard cross section has the form 



Ao-h^d(3^'<5^'Kact)GI 




1 +2(l-x) 



(4.33) 



Hence, 



/' dxAa^^^U^, Q\ /xLt)Gi = 0, (4.34) 
Jo 

and the gluonic contribution to T'l vanishes. This is so because the axial anomaly character- 
ized by the ^(1 — x) term is shifted from the hard cross section for photon-gluon scattering 
in the chiral-invariant factorization scheme to the quark spin density in the gauge-invariant 



scheme. It was first observed and strongly advocated by Bodwin and Qiu |[TT| that the above 
conclusion is actually quite general: The hard gluonic contribution to the first moment of 
Qi vanishes as long as the ultraviolet regulator for the spin-dependent quark distributions 
respects gauge invariance, Lorentz invariance, and the analytic structure of the unregulated 
distributions. Hence, the OPE result (4.24) for T\ is a general consequence of the gauge- 
invariant factorization scheme. 

We wish to stress that the quark spin density Ag'^(x) measures the polarized sea-quark 
distribution in a helicity + gluon rather than in a polarized proton. Consequently, Aq^{x) 
must convolute with AG{x) in order to be identified as the sea-quark spin distribution in a 



proton 58 



Agf (a:, /iLt) - A^f (^, /"fact) = -^(1 - x) ® AG(x, /xLt)- (4.35) 

TT 

Since the valence quark spin distribution Aq^{x) = Aq{x) — Aqs{x) is /c_L-factorization inde- 
pendent, it follows that p9[ 



Aqciix, /i^t) - Agci(a;, /iLt) = - — {1 - x) ® AG{x, /i^t), (4.36) 

TT 



which leads to 



AgGi(g') - Agci(g') = -7^ a,iQ')AGiQ'), (4.37) 
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where we have set = Q^. Eqs. (4.33) and (4.36) provide the necessary relations between 
the gauge-invariant and chiral-invariant factorization schemes. The reader may recognize 
that (4.37) is precisely the relation (3.30) obtained in the improved parton model. 



4.4 Sea-quark or anomalous gluon interpretation for F 



1 



We have seen that there are two different popular explanations for the data of Fi. In the 
sea-quark interpretation, the smallness of the fraction of the proton spin carried by the 
quarks AS = AS^, + AS^ ^ 0.30 is ascribed to the negative sea polarization which partly 
compensates the valence-quark spin component. By contrast, a large and negative sea-quark 
polarization is not demanded in the anomalous-gluon interpretation that the discrepancy 
between experiment and the Ellis- Jaffe sum rule for Ti is accounted for by anomalous gluon 
contributions. The issue of the contradicting statements about the gluonic contributions to 
the first moment of gi{x) between the improved parton model and the OPE analysis has 
been under hot debate over the past years. Naturally we would like to ask : Are these two 
seemingly different explanations equivalent ? If not, then which scheme is more justified and 
sounding ? 

In spite of much controversy on the aforementioned issue, this dispute was actually re- 
solved several years ago |]lTl. The key point is that a different interpretation for cor- 
responds to a different /^^-factorization definition for the quark spin density and the hard 
photon-gluon cross section. The choice of the "ultraviolet" cutoff for soft contributions 
specifies the factorization convention. It is clear from (3.1), (4.33) and (4.36) that to NLO 



9i{x,Q') = AgGi(x,Q2) + ^A/,(x)®AgGi(x,g^) + A<^,,(x)Gi®AG(x,Q2^ 



2 , 



Agci(x, Q') + T^Af,{x) ® Aqciix, Q') + A<^,d(a^)ci ® AG(x, 



(4.38) 



where we have set /il^^.^ = so that the InQ^/nf^^^ terms in A/g(x) and Acrh^d(x) vanish. 
As will be discussed in Sec. 7.1, the evolution of gi{x,Q^) in (4.38) is governed by the 
parton spin distributions. Therefore, the polarized structure function gi{x) is shown to be 
independent of the choice of the factorization convention up to the next-to-leading order of 
as, as it should be. This is so because a change of the factorization scheme merely shifts 
the axial-anomaly contribution between Aq{x) and A(t'^'^{x) in such a way that the physical 
proton-gluon cross section remains unchanged [cf. Eq.(3.15)]. It follows from (4.38) that 



J^g,{x,Q')dx = -(^l-^)5]AgGi(Q' 



^ (l - V) ^ (Agci(g^) - ^AG(g^)) . (4.39) 
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Hence, the size of the hard-gluonic contribution to Ti is purely a matter of the factorization 
convention chosen in defining Aq{x) and Aa'^^{x). This important observation on the kj_- 
factorization dependence of the anomalous gluonic contribution to the first moment of gi (x) 



was first made by Bodwin and Qiu (see also Manohar |T^, Carlitz and Manohar 
Bass and Thomas [jl^, Steffens and Thomas @). 



Thus far we have only considered two extremes of the fc^-factorization schemes: the 
chiral-invariant scheme in which the ultraviolet regulator respects chiral symmetry, and the 
gauge-invariant scheme in which gauge symmetry is respected but chiral symmetry is broken 
by the cutoff. Nevertheless, it is also possible to choose an intermediate factorization scheme 
which is neither gauge nor chiral invariant, so in general Agci = Ag'— AG for an arbitrary 
A (A = and A = 1 corresponding to gauge- and chiral-invariant schemes, respectively) [|61|] . 
Experimentally measured quantities do not depend on the value of A. 

Although the issue of whether or not gluons contribute to Fi was resolved six years 



ago [jTT], [T2|, the fact that the interpretation of Fi is still under dispute even today and 
that some recent articles and reviews are still biased towards or against one of the two 
popular implications of the measured Fi is considerably unfortunate and annoying. As 
mentioned in Sec. 4.1, the anomalous gluon interpretation has been deemed to be plausible 
and more favored than the sea-quark one by many practitioners in the field over the past 
years. However, these two explanations are on the same footing and all the known criticisms 
to the gauge-invariant factorization scheme and the sea-quark interpretation of Fi are in 
vain. Here we name a few: 



It has been often claimed |4^, ^ that soft contributions are partly included in 
A(Th^j(x)Gi rather than being factorized into parton spin densities because, apart 
from the soft-cutoff term, Ao"^^j(x)gi [see Eq.(4.33)] has exactly the same expression 
as (3.11) or (3.12). Therefore, the last term proportional to 2(1 — a;) arises from the soft 
region kj_ ~ << Aqcd, and hence it should be absorbed into the polarized quark 
distribution. This makes the gauge-invariant scheme pathological and inappropriate. 
However, this argument is fallacious. It is true that the 2(1 — x) term in (3.11) or 
(3.12) drops out in Acrjj'^^(x)ci because it stems from the soft kj_ region, but it emerges 
again in the gauge-invariant scheme due to the axial anomaly being subtracted from 
^'^haxd(^)ci [see Eqs. (4.31-4.33)] and this time reappears in the hard region kj_ ~ fJ^fact- 
As a result, the hard photon-gluon cross section given by (4.33) is genuinely hard ! 

A sea-quark interpretation of Ff with As = —0.10 at = 10 GeV^ has been criticized 
on the ground that a bound |As| < 0.052tQ;Q52 [^ can be derived based on the infor- 



mation of the behavior of s{x) measured in deep inelastic neutrino experiments and on 
the positivity constraint that |As(a;)| < s{x). However, this claim is quite controver- 
sial [^ and not trustworthy. Indeed, one can always find a polarized strange quark 



distribution with As ~ —0.10 which satisfies positivity and experimental constraints 
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43[| . Moreover, a sea polarization of this order is also confirmed by lattice calculations 



m, 66 



By now, we wish to have convinced the reader that it does not make sense to keep disput- 
ing which factorization prescription is correct or which interpretation is superior as they are 
equivalent. Once a set of AqGi{x), AG{x), AcrJ[^j(x)Gi or of Aqci{x), AG{x), Aal^^^{x)ci 
is chosen, one has to stick to the same scheme in all processes. 

It is worth emphasizing at this point that the equivalence of the sea-quark and anomalous- 
gluon interpretations is only applied to the first moment of gi{x), but not to gi{x) itself. 
Suppose at a certain = Ql, the data of gi{x) are reproduced either by assuming Aqs{x) ^ 
but AG{x) = in the sense of the sea-quark interpretation, or by having AG{x) ^ but 
^<ls{x) = in the sense of the anomalous gluon interpretation. It is clear that these two 
explanations are no longer equivalent at > as Aqs{x,Q'^) and AG{x,Q^) evolve 
differently. An equivalence of the first moment of gi{x) does not imply the same results for 
the higher moments of gi{x). From (4.38) it is evident that in spite of a vanishing gluonic 
contribution to Ti in the gauge-invariant scheme, it by no means implies that AG vanishes 
in a polarized proton. 

So far we have focused on the perturbative part of the axial anomaly. The perturbative 
QCD results (4.35)-(4.37) indicate that the difference Aqf^ — Ag^ is induced perturbatively 
from hard gluons via the anomaly mechanism and its sign is predicted to be negative. By 
contrast, AgP(x) can be regarded as an intrinsic sea-quark spin density produced nonper- 
turbatively. As we have emphasized in passing (see Sec. 3.1), the sea-quark helicity Agf^ 
for massless quarks cannot be generated perturbatively from hard gluons due to helicity 
conservation. The question is what is the underlying mechanism for producing an intrinsic 
negative helicity for sea quarks ? Does it have something to do with the nonperturbative as- 
pect of the axial anomaly ? The well-known solution to the Ua{^) problem in QCD involves 
two important ingredients: the QCD anomaly and the QCD vacuum with a nontrivial topo- 
logical structure, namely the ^-vacuum constructed from instantons which are nonperturba- 
tive gluon configurations. Since the instanton-induced interactions can fiip quark helicity, 
in analog to the baryon-number nonconservation induced by the 't Hooft mechanism, the 
quark- antiquark pair created from the QCD vacuum via instantons can have a net helicity. It 
has been suggested that this mechanism of quark helicity nonconservation provides a natural 
and nonperturbative way of generating negative sea-quark polarization ^ |1^]. 

There are two extreme cases for the sea-quark spin component: In one case, Aqf^lQ"^) = 
so that Aqf^ arises exclusively from the perturbative anomaly mechanism. As a result, 
AG(g2) is equal to -(27r/a,)Agfi(Q2) [cf. Eq.(4.37)] and is of order 2.5 at = 10 GeVl 
In the other extreme case, Aqf^lQ"^) = AgP((5^) so that the sea-quark polarization is 
exclusively of nonperturbative nature and AG = 0, as advocated, for example, in the chiral 
soliton model [^, The realistic case should be somewhere between these two extreme 
cases. 
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In short, the sea-quark polarization Aqf^ consists of two components: the intrinsic non- 
perturbative part Ag^ induced from the QCD vacuum via instantons and the perturbative 
part i.e., Aqf^ — Ag^ generated from the anomaly mechanism. The lattice calculation (see 
Sec. 6.1) indicates that the sea polarization is almost independent of light quark flavors and 
this suggests that it is indeed the perturbative and nonperturbative parts of the gluonic 
anomaly that account for the bulk of the negative spin component of sea quarks. 



4.5 Operator definitions for Aq and AG 

The quark spin component in the nucleon can be expressed as a matrix element of a local 
and gauge-invariant operator in the gauge-invariant fc^-factorization scheme. Since in the 
parton model Aq given by (2.15) is defined in the infinite momentum frame, we first consider 
such a frame where the nucleon is moving in the z direction with momentum = p^o oo 
and helicity +|, so that 

(Poo, T IgTS^boo, T) = AqGi. (4.40) 

This is equivalent to working in the light-front coordinate in the laboratory frame 

(p, s|g7+75g|p, s) = s+Agci, (4.41) 

where "+" is a good component in the light-front quantization formulation. It should be 
stressed that Aq is not equal to the net spin vector sum / d^plq^ip) — q^ip)] in the proton rest 
frame in the equal-time quantization formulation, where q^^^\p) is the probability of finding 
a quark flavor q in the proton rest frame with momentum p^ and spin parallel (antiparallel) 
to the proton spin Technically, the helicity and spin components of the proton are 



related to each other via the so-called Melosh transformation. The quark spin Agci is gauge 
invariant but it evolves with since the flavor- singlet axial- vector current J5 = Ql^l^Q 
has an anomalous dimension first appearing at the two- loop level The dependence 
of AEgi(Q^) will be discussed in Sec. 4.6. The evaluation of the nucleon matrix element of 
involves connected and disconnected insertions (see Fig. 4), which are related to valence 
quark and vacuum (i.e., sea quark) polarizations, respectively, and are separately gauge 
invariant. Thus we can make the identification: 

{p,s\4\p, S) = {p, S\j+\p, s)eon + {p,s\4\p, s)dis = + A#)s+. (4.42) 

q 

Interestingly, lattice QCD calculations of Aq^^ and Aqf^ became available very recently 



|5B|] . It is found that Aus = Adg = As = —0.12 ± 0.01 from the disconnected insertion 



65| . This empirical SU(3)-flavor symmetry implies that the sea-quark polarization in the 
gauge-invariant scheme is indeed predominately generated by the axial anomaly. Recall that 
sea contributions in the unpolarized case are far from being SU(3) symmetric: d> u> s. 
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connected disconnected 



Figure 4: Connected and disconnected insertions arising from the flavor-singlet axial- vector current. 

In the chiral-invariant factorization scheme one is expected to have 

(p, s\Jt\v. s) = Y: (Agci - ^Ag) s+ (4.43) 

by virtue of (4.37). The question is that can one define Agci and AG separately in terms of 
a gauge-invariant local operator ? For this purpose we write 

= - HfK" + UfK^" = + UfK'', (4.44) 

with the Chern-Simons current 

K^ = ^ e^'^P'^Alid.Al - IgUcA^Al) (4.45) 

and eoi23 = 1- Since is made of gluon fields only and Jl^ is conserved in the chiral limit 
so that {G\J^\G) = 0, it is tempted to make the identification: 

{p,s\J^\p,s) = AScis+, 

{p,s\K+\p,s) = --^AGs+. (4.46) 

ZTT 

It was originally claimed in that although the topological operator is gauge variant, 
its diagonal matrix element is nevertheless gauge invariant. The argument goes as follows. 
Consider the matrix element of K'^: 

{N{p')\K'^\N{p)) = u{p')[Gi{q')ri, + G2{q')q^lMp)- (4-47) 

Since dfj,K^ = (as/2n)GG is gauge invariant, so is the expression 2m]yGi{q'^) + q^G2{q^)- 
Consequently, the absence of a Goldstone pole coupled to implies that Gi{q^ = 0) 
and hence the matrix element of in the forward direction becomes gauge invariant. 
Another argument relies on the observation that under the gauge transformation, K^^ — >■ 
K^^ + e'^''^^dy{- ■ ■). So the gauge-dependent term can be expressed as a four-derivative and 
thus does not contribute to the diagonal matrix element of K^^. However, both above- 
mentioned arguments are erroneous for the reason of the QCD U(l) problem. In order to 
solve the U(l) problem, the SU(3)-singlet rjQ field must acquire a mass even in the chiral 
limit (see e.g., ||7l|): 

~ -Ihn^ j <fxe''^%f)\TK^{x)KM\Q). (4.48) 
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This demands a ghost pole coupled to K^. Hence, ^2(9^)9^ does not vanish in the limit 
— s> 0. Also, under the "large" gauge transformation, 

K^" ^K^ + e^'^"^9^(- • ■) + — e^'^"^Tr[(f/t9^f/)(f/t9„f/)(f/t5^f/)]. (4.49) 

It is generally believed that a solution to the U(l) problem needs two crucial ingredients: the 
axial anomaly and the instanton. The gauge transformation U (x) must be nontrivial if the 
instanton or the topological structure of the vacuum exists. It follows from (4.49) that the 
forward matrix element of K'^ is not gauge invariant under the "large" gauge transformation. 



(For an explicit example in the framework of the Schwinger model, see |]72| .) 

Since the twist-2, spin-2 gluonic operator GG (= G^^G"^'^) is gauge invariant, it has been 
proposed ||73[ to utihze the divergence equation 

= Y: 2m^qz^,q + ^GG (4.50) 

to define gauge-invariant quark and gluon spin components: 

1 _ 11- 
ASci=r {p,s\J2'^mgqi^5q\p,s), AG = {p, s\-GG\p, s) . (4.51) 

ZTflj^ g Z7Tl]\f Z 

However, this local operator definition immediately encounters several insurmountable dif- 
ficulties: for example, (i) the total light quark spin in a nucleon vanishes in the zero light 
quark mass limit, and (ii) AG and As thus defined exhibit a large isospin violation, namely 
the gluon and sea-quark spin contents of the neutron are different from that of the proton: 
Asn 7^ Asp and AG„ 7^ AGp (an explicit calculation shows AGp < 0, while AG„ > 0) 
7^ . We conclude that there is no (spin-1 or spin-2) twist-2 gauge-invariant local operator 



definition for AG and Agci- This is consistent with the OPE statement that there is only 
one twist-2 local gauge-invariant operator contributing to Fi. 

It turns out that although K^j, is not gauge invariant, its nucleon matrix element can be 
related to AG defined below in (4.56) by choosing a specific gauge and coordinate. Spin and 
orbital angular momenta in QCD are governed by a rank-3 tensor M^*^^ |5f 



+ (G'^M'^ + G'^M^) - ^G2(xV^-xV), (4.52) 

where the color indices are implicit. The fourth term in (4.52) is relevant to the gluon spin 
and the generator of gluon spin rotations has the formQ 

M°*^ (spin) = X (4.53) 



^^Note that the generators of gluon spin and orbital rotations corresponding to the respective fourth and 
third terms in (4.52), were originally incorrectly identified in |Q with Ax E and E^x x W)A\ respectively. 
However, the gluon's total angular momentum operator Jq ^ x x {E x B) given in js^ is correct. 
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However, the gluon spin and orbital terms in Mq"^ are separately gauge variant and hence 
a choice of gauge fixing is necessary. It appears that in the infinite momentum frame and in 
the temporal axial gauge A'^ = 0, the operator E x A measures the gluon spin, that is p6 



(Poo, T I X 1)' \p^, T)^o=o = AG. (4.54) 

It is easy to check that the Chern-Simons current in temporal axial gauge is proportional 
to [E X a) . We could also define the same AG in the laboratory frame using the light-front 



coordinate to obtain |[75[| 



{p,s\M+^\spm)\p, s)a+=o = (p, s| (i X i)' + 1^ ■ B^\p, s)a+=o = s+AG, (4.55) 

with Bi = ^eijkG^^, by noting that the gauge condition A" = in the infinite momentum 
frame is modified to the light-front gauge A'^ = in the light-front coordinate. Therefore, 
in light-front gauge [|, [75[ 



{p,s\K^\p,s)a+=o = -7T^AGs+. (4.56) 

ZTT 

This is the local operator definition for the gluon spin component. Consequently, we also 
have 

{p,s\4\p,s)a+=o = AScis+. (4.57) 

We see that (4.46) is valid in the light-front coordinate and in light-front gauge. 

The gluon spin AG (and likewise for ASci) also can be recast as a nucleon matrix element 
of a string-like nonlocal gauge-invariant operatoi ||76[. Of course, this nonlocal operator will 
be reduced to the local operator or Mft'^^(spin) in light-front gauge. Moreover, it is also 
possible to have operator representations for AG(a;) and Aq{x). The interested reader is 
referred to |^ |Tl| . 



From (4.57) it is clear that ASci does not evolve as the current J5 is conserved in the 
chiral limit. In the improved parton- model picture discussed in Sec. Ill, this is so because 
the ultraviolet cutoff for Aqci{x) is chiral invariant. Hence it is consistent with the naive 
intuition that the quark spin is not affected by gluon emission. Applying (4.56) and (4.57) 
to the axial-current matrix element leads to 

(p, S| J5+IP, S) = {p, s\J^\p, s)con + (P, s\J^\p, S)dis + (p, s\nfK~^\p, S)dis 

U^q? + Aq^' - ^AG)s\ (4.58) 

q ^TT 

where use of Aq^^ = Aq^^ has been made. This is in agreement with (4.43), as it should be. 
4.6 Anomalous dimensions of AE and AF 

It is pointed out in Sec. 3 that in the improved parton model there is an anomalous gluonic 
contribution to the first moment of gi{x) even in the asymptotic limit. This can be seen by 
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solving the spin- dependent Altarelli-Parisi equation (3.24). However, it can be also under- 
stood in the OPE by considering the anomalous dimension of the Chern-Simons current K^. 
The QCD evolution equation for Jg and is given by 

(4.59) 




where t = hiQ / A^^j^, and 7jj are anomalous dimensions: 

^ = ^(0) + ^^(1) + . . . . (4.60) 

/TT 

Obviously, 712 = due to the absence of Jl^ and mixing (the latter being gauge variant). 
Also, 722 = because d^K'^ ~ GG and GG does not get renormalized. Moreover, the fact 
that the Adler-Bardeen relation d^Ji^ = njdfj^K^ must be true at any renormalization scale 
fi implies that 711 = ^^721. Next consider the evolution equation -^K^ = 1^721^5* and take 
quark matrix elements. Since is of order a^, it is evident that 721 is also of order a^. As 
a result, (4.59) reduces to 




2n J \ 7(1) 




(4.61) 



Therefore, the anomalous dimension of J5 starts at the 2-loop level. The observation in 
Sec. 3.2 that a^AG is conserved to the leading-order QCD evolution is now ascribed to the 
fact that the anomalous dimension of K'^ starts at the order of and that agAG has the 
same anomalous dimension as that of since it is related to the nucleon matrix element 
of K+ via (4.56). 

Now 7^^-' can be calculated at the 2-loop level (i.e., 711) with J5 or at the 1-loop level 
(i.e., 721) with K^^. A direct calculation of 711 by Kodaira et al. gives '-y'^i = —'2nf, 
while 72}'' is computed in to be —2. Hence the relation 711 = 'n,/72i is indeed obeyed. A 
solution of the renormalization group equation 



yields 



.| + /3^+7) Jn«,9) = (4.62) 



. e.p .r(0) ^ (1 . ^M^y-) Jm (4.63, 



with Pq = (33 — 2nj)/3 and as{Q^) = 47r/(/5o In Q^/Aq^jj)). Hence, the total quark spin ASqi 
defined in the gauge- invariant fc^-factorization scheme begins evolution with at order a^. 
Since the anomalous dimension 7 is negative, AEgi((5^) decreases with Q^. 

From various operator definitions for ASqi, ASci and AF = (as/27r)AG' given in 
Sec. 4.4, it is easily shown from (4.61) that (see also ||80(| ) 

V / A v„, \ 

(4.64) 
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(4.65) 



in the gauge-invariant scheme, and 

dt\ AT ) \27rJ \2 -2n/ 

in the chiral-invariant scheme. It is evident that ASci is conserved. For parton spin densities 
Aq{x, Q"^) and AG(x, Q^), the anomalous dimensions are related to spin-dependent splitting 
functions, which we will discuss in Sec. 6.2. 



4.7 A brief summary 

It is useful to summarize what we have learned from Sees. 3 and 4. Depending on how 
we factorize the photon-gluon cross section into hard and soft parts and how we specify 
the ultraviolet cutoff on the spin-dependent quark distributions, we have considered two 
extremes of A;_L-factorization schemes. 

In the chiral-invariant factorization scheme, the ultraviolet regulator respects chiral sym- 
metry and gauge invariance but not the axial anomaly. Consequently, Aqci does not evolve 
with Q"^ and is close to the conventional parton-model intuition. There is an anomalous 
gluonic contribution to the first moment of gi{x) due to the gluonic anomaly resided in the 
box diagram of photon-gluon scattering at A;^ = [(1 — x)/Ax\Q'^ with x — > 0. Although Aqci 
cannot be written as a nucleon matrix element of a local gauge-invariant operator, a gauge- 
variant local operator definition for Aqci does exist [cf. (4.57)] in the light-front coordinate 
and in the light-front gauge — (or in the infinite momentum frame and in temporal 
axial gauge) . Since sea polarization cannot be perturbatively produced from hard gluons due 
to helicity conservation, it is expected to be small. In the extreme case that Asci — 0, AG 
is of order 2.5 at = lOGeV^, and it leads to the so-called anomalous gluon interpretation 
of Ti. 

Contrary to the above scheme, the ultraviolet cutoff in the gauge-invariant scheme satis- 
fies gauge symmetry and the axial anomaly but breaks chiral symmetry. As a result, Agoi 
is gauge invariant but dependent. Hard gluons do not contribute to Fi because the axial 
anomaly is shifted from the hard photon-gluon cross section to the spin-dependent quark 
distribution. Of course, this does not imply a vanishing AG. By contrast, the gluon spin 
component could be large enough to perturbatively generate a sizeable negative sea polariza- 
tion via the anomaly mechanism. Indeed, AG is /cx-factorization independent, and it does 
not make sense to say that AG is small in one scheme and large in the other scheme. For 
a given AG(,x), Agci and Agci are related via (4.36), which is a rigorous consequence of 
perturbative QCD. We have explicitly shown that gi{x) (not just Fi) is independent of the 
factorization prescription up to NLO. 

In order to produce sea-quark polarization for massless quarks, there are two mechanisms 
allowing for chiral-symmetry breaking and quark helicity flip: the nonperturbative way via 
instanton-induced interactions, and the perturbative way through the anomaly mechanism. 
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The empirical lattice observation of SU(3)-flavor symmetry for spin components of sea quarks 
(Sec. 6.1) suggests that it is indeed the perturbative and nonperturbative parts of the axial 
anomaly, which are independent of light quark masses, that account for the bulk of sea 
polarization. 

Although the choice of AqQi{x), Acr[|'^<j(x)Gi or Aqci{x), Aa'^^j.^{x)ci is on the same 
footing, in practice it appears that the use of AqQi{x) is more convenient than Aqci{x). 
First of all, Agci corresponds to a nucleon matrix element of a local and gauge-invariant 
operator, and its calculation in lattice QCD became available recently. For Agci, one has 
to compute the matrix element of in light-front gauge, which will require sizeable lattice 
gauge configurations. Second, NLO polarized splitting functions have been determined very 
recently in the gauge-invariant scheme, and it is straightforward to study the evolution of 
AqGi{x,Q'^) through AP evolution equations. 

5 U(l) Goldberger-Treiman Relation and Its Connec- 
tion to the Proton Spin 

5.1 Two-component U(l) Goldberger-Treiman relation 

In the gauge-invariant and chiral-invariant factorization schemes the flavor- singlet axial cou- 
pling has the expression 

9a = ASgi (5.1) 

= ASci-^AG. (5.2) 
Zn 

The smallness of the observed g\ is attributed either to the negative sea polarization or to 
the anomalous gluonic contribution. However, the question of what is its magnitude still 
remains unanswered. The well-known isotriplet Goldberger-Treiman (GT) relation 

with = 132 MeV, indicates that the coupling g\ is fixed in terms of the strong coupling 
constant g^^j^j^- It is natural to generalize this relation to the Ua{^) case to see if we can 
learn something about the magnitude of g\. 

Many discussions on the isosinglet GT relation around the period of 1989-1992 [^, ^ 
^ ^ |83i |88| , |89| were mainly motivated by the desire of trying to understand why the 



axial charge g\ inferred from the EMC experiment is so small, g\{Q) = 0.14 ± 0.17 at 
= 10.7 GeV^ (pre-1993). (The q^ of the form factor should not be confused with the 
momentum transfer occurred in deep inelastic scattering.) At first sight, the U(l) GT 
relation seems not to be in the right ballpark as the naive SU(6) quark-model's prediction 
(yf(°) = (\/6/5)g^j^j^ yields a too large value of (7^(0) = 0.80 . Fortunately, in QCD the ghost 
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field G = d^K^, which is necessary for solving the ?7a(1) problem, is allowed to have a direct 
f/yi(l)-invariant interaction with the nucleon. This together with the mixing of d^K^ with 
the ?7o implies that the net "physical" rjo — N coupling g^^j^^ is composed of the bare coupling 
9^Jnn ghost coupling g^j^j^. As a consequence, a possible cancellation between g^^j^j^ 

and g^f^^ terms will render g^ smaller. However, this two-component expression for the axial 
charge is not free of ambiguity. For example, (^^^^ is sometimes assumed to be the coupling 
between the glueball and the nucleon in the literature. Moreover, unlike the couplings g\ 
and g\, a prediction for g^ is lost. 

Since the earlier parton-model analysis of polarized deep inelastic scattering seems to 
indicate a decomposition of g\ in terms of the quark and gluon spin components 0, §], 
this has motivated many authors to identify the term {\/?>f\/2m]^)g^^j^j^ with the total quark 
spin AE in a proton, and the other term with the anomalous gluon contribution. However, 
this identification holds only in the chiral-invariant scheme. We will address this problem 
below. 

One important thing we have learned from the derivation of the isotriplet Goldberger- 
Treiman (GT) relation (5.3) is that this relation holds irrespective of the light quark masses. 
For 7^ 0, it is derived through the use of PCAC; while in the chiral limit, g\{q^) is related 
to the form factor f\{q^)q'^i which receives a nonvanishing pion-pole contribution even in the 
— > limit. By the same token, it is tempting to contemplate that the U(l) GT relation 
should be also valid irrespective of the meson masses and the axial anomaly. This is indeed 
the case: the U(l) GT relation (5.6) given below remains totally unchanged no matter how 
one varies the anomaly and the quark masses. This salient feature was first explicitly shown 
in |]8l| , |82| . It was also pointed out in p4| that this U(l) relation is independent of the 
interaction of the ghost field d^K^ with the nucleon. The easist way of deriving the U(l) 
GT relation is thus to first work in the chiral limit. Defining the form factors 

{N{p')\Jl\N{p)) = uip)[9Wh,7^ + fAiQ')Q,lMp). (5.4) 

we obtain 

2m^g'^{0) = {Nld'^j'^lN) = 3{N\d-K\N). (5.5) 
Assuming the r/o pole dominance for d-K, namely d-K = -^i^IiqI-kVO) where the r/o mass 
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arises entirely from the axial anomaly, we are led to the isosinglet GT relation]^ 



-9. 



(0) 



(5.6) 



2m7v "o"^^ 
rjo and the 

When the quark masses are turned on, chiral symmetry is explicitly broken but the GT 



with bare direct coupling between rjo and the nucleon. 



relation in terms of the r]Q remains intact, as shown in pi], R^. Nevertheless, the r]Q is no 



longer a physical meson, and it is related to the mass eigenstates via 
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where ^1, 62 and 63 are the mixing angles of tt" — 77, are given in 
values 



88 with the numerical 



-0.016, ^0 = 0.0085, 03 



48.5°. 
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In Eq.(5.7) only terms linear in small angles 61 and 62 are retained. Consequently, the 
complete GT relations in terms of physical coupling constants read 
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(5.10) 
(5.11) 

(5.12) 



^•^It is argued in jSSj g5| that the U(l) GT relation (5.6) holds only when the 7/0 is a massless Goldstone 
boson obtained in the large-A''c or OZI limit. In general the decay constant /^^ can be related to the 
topological susceptibility x'(0) of the QCD vacuum so that the U(l) GT relation reads 



2m^..90 (0) = 6./^ g^^J^ 



In the OZI limit, yx^O) = /7r/(2V3)- The smallness of the observed can be attributed either to the 
anomalously small value of the first moment of QCD topological susceptibility ^ (for an estimate of 
1/Nc corrections to x'(0), see [^) or to the suppression of the coupling gfj^J^^j^. The smallness of g\ in the 
former case is a generic QCD effect related to the anomaly and is independent of the target |Q, whereas it 
can be quite target dependent in the latter case. 

^''For the axial charge the authors of |Q obtained a result of the form (see Eq.(24) of the second 
reference of Q) 



2m N V cos 



^mv'9QNN ) - tan 6*3 ± \l ^^9^2 - Oi tan^a) 



(5.9) 



and claimed that in the limit of 6*1, 62 ^ but 9:^=^=0, it will reproduce the result of Veneziano Q] only if 
the first-order correction from ^3 (i.e., the g\03 term) is neglected. However, using Eqs.(5.11-12) and (5.16) 
one can show that (5.9) is nothing but {\/3f^/2mN)g^^J^^, as it should be. 
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where the first sign of ± or =p is for the proton and the second sign for the neutron, and the 
ellipsis in the GT relation for g'^ is related to the ghost coupling, as shown below. Since the 
mixing angles 6i and 62 are very small, it is evident that isospin violation in (5.10-5.12) is 
unobservably small. 

As we have accentuated before, the isosinglet GT relation in terms of the tjq remains 
unchanged no matter how one varies the quark masses and the axial anomaly. (A smooth 
extrapolation of the strong coupling constant from on-shell to = is understood.) 
However, the t]q field is subject to a different interpretation in each different case. For 
example, when the anomaly is turned off, the mass of rjQ is the same as the pion (for /^^ = /^). 
When both quark masses and anomaly are switched off, the rjo becomes a Goldstone boson, 
and the axial charge at = receives its contribution from the 7]q pole. 

When the SU(6) quark model is applied to the coupling g^J^^^, it is evident that the 
predicted g\ = 0.80 via the GT relation is too large. This difficulty could be resolved by 
the observation that a priori the ghost field G = d-K is allowed in QCD to have a direct 
coupling with the nucleon 

£ = ... + 3_^d^GTT{Nj,j,N) + ^{d-K)r^o + ■■■, (5.13) 

so that 

d-K = -^mlUvo + lga^^mlUd>'TT{N^,^,N). (5.14) 

However, the matrix element {N\d-K\N) remains unchanged because of the presence of the 
d-K — rjQ mixing, as schematically shown in Fig. 5 : 

{N\d-K\N) = l=f^g(^J^^-^mlUg^^^ + ^mlUg^^^ 

We see that although it is still the bare coupling g^J^^j^ that relates to the axial charge g% 
the "physical" rjo — N couphng is modified to (see Fig. 5) 

9r,oNN = S'^qLjV + '^'^IoJ'-^SgNN^ (5.16) 

where the second term arises from the r]o — d-K mixing. As a consequence, the quark model 
should be applied to q rather than to q^^^ , and we are led to 

(0) = ^^(^.oiViV - -^^lof-9GNN)- (5-17) 

This two-component expression for the U(l) GT relation was first put forward by Shore and 
Veneziano . 
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Figure 5: Contributions to the matrix element {N\d-K\N) from (1) the rjQ pole dominance, (2) a 
direct coupling of the ghost field with the nucleon, and (3) the d-K — rjo mixing. 



It has been proposed that the smallness of Qj^ may be explained by considering the pole 
contributions to d-K from higher single particle states X above the rjo, so that the isosinglet 
GT relation has the form (see e.g., Chao et al. ||8^, Ji [^, Bartelski and Tatur ||90|1 ) 



^a(O) = ^^^fvo9,,NN + Y.fx9xNN)- (5-18) 

The state X could be the radial excitation state of tiq or a ^ glueball. (Note that the 
ghost field d-K is not a physical glueball as it can be eliminated via the equation of motion.) 
However, we will not pursue this possibility further for two reasons: (i) It is entirely unknown 
whether or not the X states contribute destructively to (7^. (ii) As we shall see later, the 
contribution from a direct interaction of the ghost field with the nucleon corresponds to 
a disconnected insertion, which is shown to be negative according to recent lattice QCD 



calculations ||6^, Therefore, the ghost-field effect is realistic, and if the contributions due 
to the states X are taken into account, one should make the following replacement 

1 2 1 2 

driQNN *' driQNN ^^^»?0 f'^dcNN ' 9xNN * 9xNN 'yj^x9 XNN (5.19) 

in Eq.(5.18), where A is the d-K~X mixing. 
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5.2 Interpretation of the U(l) Goldberger-Treiman relation 

By comparing (5.17) and (5.2), it is tempting to identify the two components of the U(l) 
GT relation as 

However, this identification is not unique and sensible because it does not hold in the 
gauge-invariant factorization definition for Aqqi. One may ask can one have a physical 
interpretation valid for both /c^-factorization schemes for the g^^^^,^ and g^,^,^ terms in the 
two-component isosinglet GT relation (5.17) ? As noted in Sec. 4.5, the evaluation of the 
hadronic flavor- singlet current involves connected and disconnected insertions (see Fig. 4) 
which are related to valence-quark and sea-quark contributions respectively and are sepa- 
rately gauge invariant. A recent lattice calculation indicates an empirical SU(3)-fiavor 
symmetric sea polarization; this implies that the disconnected insertion is dominated by the 
axial anomaly of the triangle diagram. Since the triangle contribution is proportional to 
d-K, the ghost field, it is thus quite natural to make the gauge-invariant identification: 



73/. _ 

In """" 2mAr 



connected insertion, — 7r^~'9GNN — disconnected insertion, (5.21) 



which is valid in both factorization schemes. Note that this identification is basically an 
assumption since it is possible to add and substract some part of the disconnected contrib- 
tuion in (5.21) and the resultant identification is still gauge invariant. In the gauge-invariant 
factorization scheme, the disconnected insertion, which is responsible for the smallness of 
g% should be interpreted as a screening effect for the axial charge owing to the negative sea 
polarization rather than an anomalous gluonic effect. 

Having identified the two-component U(l) GT relation (5.17) with connected and discon- 
nected insertions, we are now able to extract the physical coupling constants g^,j^^ and g^^^^ . 
This is because the connected insertion (CI) corresponds to the total "valence" quark con- 
tribution (strictly speaking, the valence-quark plus cloud-quark contributions; see Sec. 6.1) 
to the proton spin, so it is related to the quark model expectation; that is, 

= 9a{CI) = Au. + =3F-D, (5.22) 

where the last identity follows from the fact that g\ = 3F — D = Au + Ad — 2As 
Auy + Ady due to the aforementioned SU(3) symmetry for sea polarization. Unlike the 
previous identification (5.20), g^{Cl) here is not identified with the total quark spin AE. 
In the non-relativistic quark limit, F = |, D = 1, and hence Au^ + Ad^ = 1. With the 
inclusion of the relativistic effects and cloud-quark polarization (see Sec. 6.1), F and D are 
reduced to 0.459 and 0.798, respectively, and (7^(CI) is reduced to a value of 0.579 . 
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From Eqs.(5.11), (5.12) and (5.22), the GT relations for g\ and are recast to 
or — V = — fif = — (fif cost^3 + fif ,„„ smt^s), 

— D = ^fif .... = ^(fif , cos6'3 — fif sin ^3), (5.23) 

2m n""" 2mN """"^ " 

where the tiny isospin-violating effect has been neglected. Note that we have fif^^^^jv instead 
of on the second line of the above equation. Using 63 = —18.5° [see Eq.(5.8)], it follows 
from (5.23) that [|T| 

W=3-4, ^?,^^=4.7, (5.24) 

while 

^7.0--= 4-8. ^7.3iviv=3.4. (5.25) 

It is interesting to note that we have g^,^^ < g^,^^, whereas ^f^^^^^ > g^^^^j^- Phenomenologi- 
cally, the determination of g^,^j^ and g^^^^^ is rather difficult and subject to large uncertainties. 
The analysis of the A^A^ potential yields g^,^j^ = 7.3 and g^^^^ = 6.8 while the forward 
NN scattering analyzed using dispersion relations gives g^,^^^ g^r^j^ < 3.5 |9^. But these 



Finally, the ghost coupling is determined from the disconnected insertion (DI) 

-9gmn = ^7° (DI) = Aus + Ads + As ^ 3As. (5.26) 



analyses did not take into account the ghost pole contribution. An estimate of the rj' 27 
decay rate through the baryon triangle contributions yields g^,j^^ = 6.3 ± 0.4 

ing is determined from the discom 

2mN 

Using g^^{0) = AS = 0.31 ± 0.07 [see (2.27)] and (5.26) we obtain 

g^,,^55GeV-\ (5.27) 

In principle, this coupling constant can be inferred from the low-energy baryon-baryon scat- 
tering in which an additional SU(3)-singlet contact interaction arises from the ghost inter- 
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To summarize, the U(l) GT relation (5.6) in terms of the t]o remains totally unchanged 
no matter how one varies the quark masses and the axial anomaly, while its two-component 
expression (5.17) can be identified with the connected and disconnected insertions. Since 
(-\/3/^/2mAr)(yf^^j^j^ is related to the total valence quark contribution to the proton spin, we 
have determined the physical coupling constants g^,^^ and g^j^^ from the GT relations for 
9a and g^^ and found that g^,^^ = 3.4 and ^f^^^ = 4.7 . 



6 Other Theoretical Progresses 
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6.1 Lattice calculation of proton spin content 

The spin- dependent DIS experiments indicate that Am ~ 0.83, Ad ~ —0.43 and As ~ —0.10 
at = 10 GeV^ [cf. (2.26)]. We learn from Sees. 3 and 4 that the axial anomaly plays 
an essential role for the smallness of g\ or the suppression of Fi relative to the Ellis- Jaffe 
conjecture. However, many questions still remain unanswered, for example: (i) what is 
the sea polarization of the non-strange light quarks (i.e., Aug, Adg) ? (ii) what are spin 
components of valence quarks Au^, Ad^ ? Are they consistent with the expectation of quark 
models ? (iii) what is the magnitude and sign of the gluon spin component in a proton ? 
(iv) what is the orbital angular momentum content of quarks and gluons ? and (v) what are 
spin-dependent parton distributions Ag(x), AG{x) ? A truly theoretical or experimental 
progress should address some of the above-mentioned questions. Obviously, a first-principles 
calculation based on lattice QCD will, in principle, be able to provide some answers. Indeed, 
the present lattice calculation is starting to shed light on the proton spin contents. 

After the 1987 EMC experiment, there existed several attempts of computing AG and g\ 
using lattice QCD (for a nice review, see Liu Okawa and references therein). A first 
direct calculation of the quark spin content Aq was made in but without final results. 
Fortunately, two successful lattice calculations in the quenched approximation just became 



available very recently ||6^, A more ambitious program of computing the polarized 

structure functions gi{x), g2{x) and their moments is also feasible and encouraging early 
results were reported in PB| . 

What computed in ^ is the gauge-invariant quark spin component Agci defined 
by s^AqQi = {p, slq'j^'y^qlp, s) (recall that Aq has the conventional partonic interpretation 
only in the component in the light-front coordinate). An evaluation of Agci involves a 
disconnected insertion in addition to the connected insertion (see Fig. 4; the infinitely many 
possible gluon lines and additional quark loops are implicit). The sea-quark spin contribution 
comes from the disconnected insertion. It is found that 

[65]: Audis = A4is = -0.12 ±0.01, As = -0.12 ± 0.01 , 

H : Awdis = Ac^dis = -0.119 ±0.044, As = -0.109 ± 0.030 . (6.1) 



Note that the results of fS^] are gauge independent although the gauge configurations on 
the t = time slice are being fixed to the Coulomb gauge (see a discussion in [Q). It is 
evident that the disconnected contribution is independent of the sea-quark mass in the loop 
within errors. Therefore, this empirical SU(3)-fiavor invariance for sea polarization implies 
that the disconnected insertion is dominated by the axial anomaly of the triangle diagram; 
that is, it is the gluonic anomaly which accounts for the hulk of the negative sea polarization. 
This is consistent with the picture described in Sec. 4.3, namely a substantial polarization 
of sea quarks is produced from gluons via the perturbative anomaly mechanism and from 
nonperturbative effects via instantons. 



45 



It has been emphasized in [£9| that the connected insertion involves not only valence 
quarks but also cloud quarks. In the time-ordered diagrams, one class of the connected 
insertion involves an antiquark propagating backward in time between the currents and is 
defined as the "cloud" antiquark as depicted in Fig. 6. Another class involves a quark 
propagating forward in time between the currents and is defined to be the sum of valence 
and cloud quarks. Hence the quark spin distribution can be written as 

Aq{x) = Aqvix) + Aq,{x) + Ag,(x) = Aq^{x) + Ag,(x), (6.2) 

where Aqy{x) as conventionally referred to as the "valence" quark spin density is actually a 
combination of cloud and truly valence contributions, i.e., Aq^{x) = Aqv{x) + Aqc{x). The 
concept of cloud quarks, which is familiar to the nuclear-physics community, appears to be 
foreign to the particle-physics community. As shown in [^, the presence of cloud quarks 
and antiquarks is the key for understanding the origin of deviation of the Gottfried sum rule 
from experiment, namely the difference of u and d distributions in the nucleon. 





Figure 6: Time-ordered diagrams of the connected insertion involving quark and antiquark prop- 
agators between the currents. 

A very important lattice observation is made in that the SU(6) relation in the quark 
model is recovered in the valence approximation under which cloud quarks in the connected 
insertion are turned off. For example, the ratio Ra = fl'A.con/S'i is found to be | in the lattice 
calculation when the presence of cloud quarks and antiquarks is eliminated by disallowing 
quarks from propagating backward in time, while theoretically it is reduced under the valence 
approximation to (A-Ucon + A(icon)/(A'Wcon — ^c^con), which is equal to | in the relativistic 
or non-relativistic quark model. Of course, the prediction (5'^)'^^ = | in the non-relativistic 
quark model is too large compared to the experimental value (5'^)expt = 1.2573 ±0.0028 



Presumably, ((7^)'^^ is reduced by a factor of 3/4 due to relativistic effects. In other words, the 
above lattice observation implies that relativistic quark model results should be recovered in 
the valence approximation. Based on this observation, a picture for the smallness of AEqi or 
g\ emerges. In the relativistic quark model, the non-relativistic SU(6) predictions Am^^ = | 
and Afi"^^ = — | are reduced by the same factor to Auy = 1 and Ady = — i, where the 
subscript "V" denotes a genuine valence spin component. Since the quark orbital angular 
momentum is nonvanishing in the presence of quark transverse momentum in the lower 
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Table II. Axial couplings and quark spin contents of the proton from lattice calculations and from 
experiments [see (2.19) and (2.26)]. 





H 


11 


Experiment 


9'a 


0.25(12) 


0.18(10) 


0.31(7) 


9\ 


1.20(10) 


0.985(25) 


1.2573(28) 


9\ 


0.61(13) 


— 


0.579(25) 


Au 


0.79(11) 


0.638(54) 


0.83(3) 


Ad 


-0.42(11) 


-0.347(46) 


-0.43(3) 


As 


-0.12(1) 


-0.109(30) 


-0.10(3) 


F 


0.45(6) 


0.382(18) 


0.459(8) 


D 


0.75(11) 


0.607(14) 


0.798(8) 



component of the Dirac spinor, the reduction of the spin component from Am^^ + Ad^^ = 1 
to Auy + Ady = 0.75 is shifted to the orbital component of the proton spin. Assuming 
SU(3)-symmetric sea polarization, as suggested by lattice calculations, one obtains from 
(2.26) that 

Au^ = Auv + Auc ~ 0.93 , Ad^ = Ady + A4 ~ -0.33 , (6.3) 

and hence 

Auc ~ -0.07 , A4 ~ -0.08 . (6.4) 

The cloud-quark polarization Aq^ is thus negative in sign and comparable in magnitude to 
the sea polarization Aqg. Now we have 

Amgi + Ada + AgQi = Auv + Ady + Au^ + A4 + Aus + Ads + Ag . (6.5) 

0.30 0.75 -0.15 -0.30 

We conclude that the deviation o/AEgi or g\ from the relativistic quark model's value 0.75 
is ascribed to the negative cloud-quark and sea-quark polarizations. In the future, it will be 
of great importance to calculate Aqy and Aq^ directly by lattice QCD. 

The lattice results of |6^ and for Agci, Qa, F and D are presented in Table II; in 



general they agree with experiments within errors. Moments of polarized structure functions 
gi{x) and (72 (x) from the connected insertion are reported in It is found that twist- 



3 operators characterized by the matrix element ^3 [cf. Eq.(4.18)] provide the dominant 
contribution to x'^g2{x)dx. 

As for the chiral-invariant quantity ASci, it involves the matrix element of in light- 
front gauge [see Eq.(4.57)] and hence sizeable gauge configurations are needed in lattice 
calculations for AEci- Nevertheless, it is conceivable to have lattice results for AG and Agci 
soon in the near future. 
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6.2 Two-loop spin-dependent splitting functions 

The experimental data of gi{x^Q'^) taken at different x-bin correspond to different ranges 
of Q^, that is, of the data is x-bin dependent. To the zeroth order in QCD, gi simply 
reads gi{x) = \ Agj(x) without scaling violation. To the leading order (LO), it becomes 
gii^x^Q"^) = I Z^i Agj(x, Q^) with scaling violation arising from gluon bremsstrahlung and 
quark- antiquark pair creation from gluons. In other words, AG{x) enters into gi at LO only 
via the evolution governed by the LO polarized AP equation. To the next-to-leading order 
(NLO), gi{x, Q"^) is given by (3.1). At this order, gluons contribute directly to the polarized 
structure function gi. A full NLO QCD analysis of the gi data is thus not possible until the 
two-loop splitting functions APj^j\x) in the NLO evolution equation are known. Since the 



complete results for APj^j\x) are not available until very recently [^], all pre-1995 analyses 
based on the NLO expression (3.1) for gi{x,Q'^) are not complete and fully consistent. 

The dependence of parton spin densities is determined by the spin-dependent Altarelli- 
Parisi equations: 



^AgNs(x,t) = ^AP^7(x) ® AgNs(x,t), 

d ( Aqsix, t) \ _ «,(t) / Aps^(x) 2nfAP,Gix) \ ( Aqs{x, t) 



(6.6) 



with t = ln(QVA2 



dt\AG{x,t)J 271 \APGq{x) APgg{x) J \AG{x,t) J 
^ ) 

QCD''' 

Aqj,s{x) = Aqi{x) - Aqjix), Ags(a;) = J2Aqi{x), (6.7) 



and 



AP,,(x) = Ap!^\x) + ^Ap(^)(x) + ■ ■ ■ . (6.8) 



The spin-dependent anomalous dimensions are defined as 



A7I;. = AP,,(x)x"-Mx = A7jf " + gA7,y)'" + ■ ■ ■ . (6.9) 



The leading-order polarized splitting functions AP/J'^ are given by (3.26) and the correspond- 
ing anomalous dimensions for n = 1 are 

A7r = ^7jr = 0' A7a^' = 2, ^lPa' = lPo, (6.10) 
where /?o = H - 2n//3. To the NLO, AP^W and AP^^^^ were calculated in the MS scheme by 



Zijlstra and van Neerven [|100|| . However, the other two polarized splitting functions APq^ 



and ApS were not available until last year. The detailed results of AP^^\x) are given in 



51| . Here we just list the anomalous dimensions for n = 1: 
A7{^^i = 0, A7^;],^ = -3C^T^ = -2n^, 

= 0' ^7^'^ = -I {-QCl - '-^C^C, + l^T,) = 25 - Inj, (6.11) 
(i),i _ A _ If 68 2 , _ _ , 40^ _ \ 1 A_ 38 



^i'gs = X = -8 ry + ^^^^^ + y^^^^ ) = i v^'^ " y""^ 
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where Cp = |, Ca = 3, Tj = nf/2. To this order, 

1 AlnlngVAg, 
47r /3olnQVA^ Pli\nQyAl_y 



(6.12) 



Note that the MS regularization scheme is a gauge-invariant fc^-factorization scheme as 
it respects the axial anomaly in the triangle diagram. Therefore, the NLO evolution of par- 
ton spin distributions in the gauge-invariant factorization scheme is completely determined. 
Explicitly, the AP equation for the first moment of flavor- singlet parton spin densities reads 



■ASGi(t)\ ^ «^ / t(-2n;) \ /ASGi(t) 

AG{t) ) 2n [2 + l^{25 - luf) f + f A ^G{t) 

Defining AF = (Q;s/27r)AG, it is easily seen that 



(6.13) 



/ASGi(t)\ / OWASGi(t) 

dt [ AT{t) ) \2n) [2 0) [ AT{t) 



0{al), (6.14) 



which is in agreement with (4.64). A derivation of (6.14) does not need the information 
of A7^^g'^ and A'-^qq^^ however. The NLO evolution of parton spin densities has been 
studied in [|62| , |101| , |102| , |103|| . It is found that the difference between LO and NLO evolution 
for AIl(x, Q^) and AG{x, Q"^) is sizeable at small a;, a; < 5 x 10~^ (see Figs. 4 and 5 of ||103|| ). 
This feature can be understood from the x — > behavior of the splitting functions APjj(x). 
As X — 0, we find from (3.26) that 

APi;)-t APig^-l, AP§l^l AP(S^12, (6.15) 



and from [511 that 



AP^^) ^ (aCfCa - SCpTf - 6Cl) In^ a; = y (1 - nj) In^ x, 

22 

AP^(^) ^ - {2Ca + Cp) In' x = -- In' x, 



AP^al (^CaCf + 4C|) In^ x = — In^ x, (6.16) 

16 



AP^CG - ^CpTf) \n' X = y(27 - rij) \n' x. 

It is evident that for small enough x, the NLO AP//^ (x) can overcome the suppression factor 
(as/27r) and become comparable to the LO splitting functions. 



6.3 Orbital angular momentum 

We have discussed the operator definitions for Aq and AG that are accessible with experiment 
in Sec. 4.5 and their Q' evolution in Sees. 3.3 and 4.6. It is natural to see if the similar 
analysis can be generalized to the orbital angular momenta of quarks and gluons. So far 
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we have noticed two places where the orbital angular momentum plays a role. One is the 
compensation of the growth of AG with by the angular momentum of the quark-gluon 
pair (see Sec. 4.3). The other is the reduction of the total spin component AE due to the 
presence of the quark transverse momentum in the lower component of the Dirac spinor is 
traded with the quark orbital angular momentum (see Sec. 6.1). 
The generators associated with rotation invariance are 

Jf'^ = J d^xM^^"", (6.17) 

where M"'^^ is the angular momentum density given by (4.52). The first and third terms 
in (4.52) contribute to the quark and gluon orbital angular momentum, respectively. The 
angular momentum operator in QCD is related to the generators by 

Exphcitly H], 

Jl = Sl + Ll = j (fx [^t'tsV' + i^\xx idfi)], 

Jf = + L^^ = j (fx[{Ex Af -E,{xxdfAi]. (6.19) 
Except for the quark helicity operator SI, the other three operators L^, S*^, are not 



separately gauge and Lorentz invariant. Very recently, Ji [104] has obtained gauge- invariant 
expressions: 

Ll= [ dfxip\x X iDfip, Jf = f dfx [x X (^ X B)f . (6.20) 



The gluon total angular momentum Jf does not permit further gauge-invariant decompo- 
sition into spin and orbital pieces. However, in the infinite momentum frame and in the 
temporal axial gauge A^ = 0, measures the gluon spin component AG which is accessible 
experimentally [cf. Eq.(4.54)]. As a result, the nucleon matrix element of in the infinite 
momentum frame and in A" = gauge (or in the light-front coordinate and in light-front 
gauge) can be deduced from the matrix elements of and ||104|| . However, whether this 
definition of the gluon orbital angular momentum (and likewise Ll) contacts with experiment 
is still unknown. 

The evolution of the quark and gluon orbital angular momenta was first discussed by 
Ratcliffe [|105|| . Using the operators given in (6.19), Ji, Tang and Hoodbhoy [0] recently 



have derived a complete leading-log evolution equation: 



dt\LfJ 27r \ ICf J KL'^ J 27r \-ICf 




with the solutions 



= + - — ^ + f{Q') + -AE 



2 2 16 + 3n/ 2 2 16 + Sri/ ^ 
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L<:{Q^) = -AGiQ^) + l-^ + f{Q') L<:{Ql) + AGiQ'^ ^ 



2 16 + 3nf ' ■"'^ ' \ ' ' ^""^^ 2 16 + 3/2/^ 



(6.22) 



where 



'inQg/A 



2 



and AS is independent to the leading-log approximation. We see that the growth of 
AG with Q"^ is compensated by the gluon orbital angular momentum, which also increases 
like InQ^ but with opposite sign. The solution (6.22) has an interesting implication in the 
asymptotic limit ^ oo, namely 



2 ' 2 16 + 3n/ ' 

1 1 6 

J^{Q') = AG{Q') + L^{Q')^ ^Y^:^^- (6.24) 

Thus, history repeats herself: The partition of the nucleon spin between quarks and gluons 
follows the well-known partition of the nucleon momentum. Taking = 6, we see that 
J| : = 0.53 : 0.47 . If the evolution of and is very slow, which is empirically 
known to be true for the momentum sum rule that half of the proton's momentum is carried 
by gluons even at a moderate Q^, then AS ~ 0.30 at = 10 GeV^ implies that L« ~ 0.10 
at the same Q^, recalling that the quark orbital angular momentum is expected to be of 
order 0.125 in the relativistic quark model. 

Finally, it is worthy remarking that the spin sum rule 

i = 1 ASgi + (L^)gi + AG + (6.25) 

so far is defined in the gauge-invariant fc^-factorization scheme. In the chiral-invariant scheme 
we have ASci = ASqi + {n fas/ It:) AG, but gi{x) and Fi remain unchanged [cf. Eq.(4.38)]. 
Since AG and are independent of the fc^-factorization, a replacement of ASqi by ASci 
in the spin sum rule (6.25) requires that the difference ASci — ASgi = {nfas/27i)AG be 
compensated by a counterpart in the gluon orbital angular momentum; that is, 

11 -fry 

{LDci = {Llhi - ^AG. (6.26) 



This relation also can be visualized as follows (see [§2|). Suppose we first work in the 
chiral-invariant scheme and consider a gluon with +1 helicity splitting into a massless quark- 
antiquark pair. The total helicity of the gluon is entirely transferred to the orbital angular 
momentum of the pair due to helicity conservation or chiral symmetry. Now, shifting the axial 
anomaly from the hard part of the photon-gluon box diagram to the triangle diagram so that 
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a negative sea-quark polarization is produced via the anomaly mechanism [see Eqs.(4.29)- 
(4.32)]. In order to preserve the total angular momentum, this sea-quark polarization must 
be balanced by the same amount of the quark orbital angular momentum induced from the 
anomaly. It is interesting to note that when AG is of order 2.5 , one will have ASci ~ 0.60 
[cf. Eq.(3.31)] but {LDci ~ —0.05. In other words, while ASci is close to the quark-model 
value, {Ll)ci deviates more from the quark model and even becomes negative ! 



7 Polarized Parton Distribution Functions 
7.1 Prelude 

One of the main goals in the study of polarized hadron structure functions measured in 
DIS is to determine the spin-dependent valence-quark, cloud-quark, sea-quark and gluon 
distributions and to understand the spin structure of the nucleon. In spite of the recent 
remarkable progress in polarized DIS experiments, the extraction of spin-dependent parton 
distribution functions, especially for sea quarks and gluons, from the measured polarized 
hadron structure functions remains largely ambiguous and controversial. We shall see that 
a full NLO analysis of the gi{x,Q'^) data just became possible recently and it indicates 
that the sea-quark and gluon spin distributions are, to a large degree, still unconstrained 
by current experimental data. Nevertheless, we are entering the phase of having the parton 
spin densities parametrized and determined to the NLO. 

In general the polarized proton structure function gi{x, Q^) has the form [|100|| 

9i{x,Q') = \Y.^l[^^{x,^4.c^ + '^^^^[^^^^ (7.1) 

-(NS/^ ,n2 ,.2 \^ \ „ 1^ ,.2 



with 



AQg(x) = ACj^^Ux) + ^AC%{x) + ■ ■ • . (7.2) 

Note that (as/27r)ACG(x) is equal to the hard photon-gluon cross section Acr^^^(a;) in (3.1) 
and ® denotes convolution. The gluon coefficient function ACg{x) and the quark spin 
density Ag(x) depend on the fc^-factorization scheme, while the quark coefficient function 
ACq(x) depends on the regularization scheme chosen. In the MS scheme, which is also 
a gauge-invariant factorization scheme, ACg°](a;, Q^) = AC^l^q{x,Q'^) = A/^(x, Q^) [cf. 
Eq.(3.3)] and 



ACPix, Q\ /iLt)Gi = (2x - 1) (in + In ^ - l) + 2(1 - x) (7.3) 
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and the NLO AC^^^ are given in ||100|| . In the chiral-invariant scheme, the leading order 
ACg{x) is calculated to be [cf. Eq.(3.20)] 

ACiS\x, Q\ /iLt)ci = (2a: - 1) fin ^ + In - 1 V (7.4) 

The first moments of the coefficient functions are 

ACf\x)Gi = -2, £ a4°)(x)gi = 0, £ a4°)(x)ci = -1. (7.5) 

The dependence of the parton spin densities is determined by the AP equation (6.6). As 
mentioned in Sec. 6.2, at the zeroth order of a^, ACq^cix) = and APij{x) = 0. At NLO 
we still have ACg^cix) = but AP^j\x) ^ 0; that is, there is a scaling violation in gi{x, Q"^) 
but AG{x) enters indirectly via the evolution. A complete NLO analysis of gi{x,Q^) 
requires the information of AP^^\x) in addition to AC^^q{x). At this order, gluons start to 
contribute directly to the polarized structure function. For a next-to-next-to-leading order 
description we have to await three-loop results for APj^j\x), although ACjj^Q{x) have been 
calculated. 



Several remarks are in order. 



It is clear from (7.1) that with the input of parton spin distributions at = fif^, 
the evolution is governed by the logarithmic term ln((5^//Ufact) ^^e coefficient 
functions, as long as ^^(/ifact) lii(Q^//ifact) << 1- » /^facf ^^e logarithmic 

terms have to be resummed using renormalization group methods |p.00|| . For a fixed 



act 1 



/ifact ^iid '^ot deviating too much from ^f^^ty ^^e ln(Q^//i£a^J terms in ACjj^Q give 

rise to the leading-log (LL) Q"^ evolution to (?i(x, Q^), and the ln^(Q^//ifaj,J terms in 
ACg^Q determine the dependence to the next-to-leading log (NLL) approximation. 
When /ifact is set to be Q^, the ln(Q^//i£apJ terms appearing in coefficient functions are 
equal to zero and (7.1) becomes 

9ii^,Q') = lj:ellAq{x,Q') + ^4r^[AC'g{x,as)®Aqs{x,Q') 



2V n ^ ' 27r 
+ ACf (x,a,)® AgNs(a:,Q') + ACG(x,a,)® AG'(x,Q2)1|. (7.6) 



In this case, the evolution of gi{x, Q"^) is taken over by the parton spin distributions. 



Using F2{x, Q ) as a testing example, it is shown explicitly in ||100|| that the Q depen- 



dence determined by the LL (NLL) parametrization of parton densities in which LL 
(NLL) logs are resummed to all orders of perturbation theory is indeed consistent with 
the leading (next-to-leading) evolution obtained from fixed order perturbation the- 
ory (i.e., /ifact being kept fixed). In short, generally we have to solve the spin-dependent 
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AP equation (6.6) to determine the dependence of spin- dependent parton distribu- 
tions and hence the evolution of g\ via (7.6). However, for not deviating too 
much from //fact) ^ provides a good approximation to the evolution of g\{x, Q"^) 
through the ln(Q^//i£a^^) terms in coefficient functions. 

Although the contribution ACg^AG in (7.1) or (7.6) is formally of order a^, it actually 
does not vanish in the asymptotic limit due to the axial anomaly. It is thus expected 
that the NLO corrections to sea-quark and gluon spin distributions are important. 

Before the availability of the two- loop splitting functions APj^j\x), some analyses of 
gi{x,Q^) were strictly done at the leading order, namely gi{x,Q'^) = | e^^g(x, Q^) 
with the gluon-spin effects entering via the evolution (see e.g., [|106| , 0). As 



as{Q^)AG{Q^) is of order several analyses have been performed using a hybrid 
expression for gi 



gi{x,Q^) = 



^q{x, Q') + ^ ACg{x) ® AG{x, 



(7.7) 



in the chiral-invariant factorization scheme. However, the gluon coefficient function 
employed in many earlier studies is often incorrect. For example, ACq^{x) = —6{l — x) 
was used in |10^ and ACgi(x) = (2x - 1) ln[(l - x)/x] in |T09| . 



In spite of the fact that the combination Aq{x, /ifact) + (as(/ifact)/2vr)ACG® AG(x, /ifact) 
in (7.1) is fc_|_-factorization independent [see Eqs.(4.33) and (4.37)], the lack of knowl- 
edge on the splitting functions APij\x) in the chiral-invariant scheme indicates that, in 
practice, we should work entirely in the gauge-invariant factorization scheme in which 
hard gluons do not make contributions to Fi. This is further reinforced by the observa- 
tion that in the literature most of NLO parametrizations of unpolarized parton distri- 
butions, which are needed to satisfy the positivity constraints \Aq{x,Q'^)\ < qi^x^Q"^) 
and |AG(x, Q^)| < G(x, Q^), are performed in the MS scheme. 



7.2 Constraints on polarized parton distributions 

As stressed in Sec. 6.1, the quark spin density Aq{x) consists of valence-quark, cloud-quark 
and sea-quark components: Agy(x), Aqc{x) and Aqs{x). Unfortunately, there is no any 
experimental and theoretical guidehnes on the shape of the spin-dependent cloud-quark dis- 
tribution, though it is argued in Sec. 6.1 that the cloud-quark polarization is comparable 
to sea-quark polarization in sign and magnitude. Since the spin component of cloud quarks 
originates from valence quarks (there is no cloud strange quark), we will proceed by consider- 
ing the combination Aq^^x) = Aqy^x) -\- Aqd^x), which is commonly (but not appropriately) 
referred to as the "valence" quark contribution. Since the sea polarization is found to be 
SU(3)-flavor symmetric empirically in lattice calculations we will make the plausi- 

ble assumption of SU(3)-symmetric sea-quark spin components. This assumption is justified 
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since the disconnected insertions from which the sea-quark spin component originates are 
dominated by the triangle diagram and hence are independent of the hght quark masses 
in the loop.^^ Therefore, for SU(3)-symmetric sea polarization, we obtain from (2.26) that 
Au^ = 0.93, = -0.33 at = 10 GeV^ [cf. Eq.(6.3)]. As explained in Sec. 6.1, the 
deviation of the result Au^ + Ad^ = 0.60 from the relativistic quark model's prediction 
Auy + Ady = 0.75 stems from the negative cloud-quark polarization. 

The valence quark spin density at x ^ 1 is subject to a model-independent constraint. 
According to the perturbative QCD argument [ |110| ], the valence quarks at x = 1 remember 
the spin of the parent proton, i.e., Auv{x)/uv{x), Ady{x)/dv{x) ^ 1 as a; — > 1, as originally 
conjectured by Feynman ||111||. Since Arf^, is negative while Adv{x) is positive as a; — > 1, 



it means that the sign of Adv{x) flips somewhere between < x < 1 ||112|] . A model for 
valence-quark spin distributions has been proposed sometime ago by Carlitz and Kaur ||113|] . 
According to this model, Adv{x)/dy{ 1, which disagrees with what expected 

from perturbative QCD. Experimentally, it is possible to carry out a straightforward mea- 
surement of the ratio Ad^{x) / dy{x) to test various predictions by measuring the longitudinal 
spin asymmetry in the inclusive W~ production in proton-proton collisions ||114|| . This spin 
asymmetry is proportional to Adfj{x)/dj;{x) in the appropriate kinematic range (see Sec. 8). 
In terms of valence and sea spin distributions, gf can be recast to the form 




Au^{y,Q^) + -Ad^{y,Q 



9 



+ -Asiy,Q'] 



{7.t 



+ 



ACi^\x/y)AG{y,Q') 



In general, both sea quarks and gluons contribute to gi{x). Since the unpolarized sea distri- 
bution and the unpolarized gluon distribution multiplied by a^/ (67r) are small at x > 0.2, 
the positivity constraints |As(x)| < s(x) and |AG'(x)| < G{x) imply that the data of gi{x) 
at X > 0.2 should be almost accounted for by Am„(x) and Ady{x). Therefore, the shape 
of the spin-dependent valence quark densities is nicely restricted by the measured gi{x) at 
X > 0.2 together with the first-moment constraint (6.3) and the perturbative QCD require- 
menlf^ that Aqv{x) / qv{x) ^ 1 at x = 1. In order to ensure the validity of the positivity 



condition |Ag^(x)| < g^(x), we choose the NLO Martin-Roberts-Stirling MRS(A') set ||116 
parametrized in the MS scheme at = 4 GeV^ as unpolarized valence quark distributions 



m,(x,Q2 = 4GeV2) = 2.26x-°-^^^(l 



x?-''Hl 



0.54v/x + 4.65x), 

^■^ It was first noticed in |101| that even if sea-quark polarization is SU(3) symmetric at, say 



Qt the 

SU(3)-flavor symmetry will be broken at > Qq due to a nonvanishing NLO Ajqq'^ in (6.11). However, 
the degree of SU(3) breaking is so small that we can neglect it. 

^^It was assumed in |5|, 115] that Auy{x) = a{x)uy{x), Adv{x) = f3{x)dv{x) with a{x), p{x) — > 1 as 
1 and a{x), P{x) — ^ as a; ^ 0. However, the constraint at a; = is not a consequence of QCD. In the 

0.136 at a; = 0. As a result, |A(ji,(a;)| 



a; ^ i ana a[x), p[x) — ^ as a; 

present work we find that Auy{x)/uy{x) = 0.41 and Ady{x)/dy{x) 
is usually larger than |As(x)| even at very small x. 
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4(x,Q^ = 4GeV2) = 0.279x-°-^^^(l-x)^-^''(l + 6.80yx + 1.93a;). (7.9) 

Accordingly, we must employ the same MS scheme for polarized parton distributions in order 
to apply the positivity constraint. For the spin-dependent valence distributions we assume 
that they have the form 



Aq^{x) = - xf{a + + cx + dx^-^), (7.10) 

with a and (3 given by Eq.(7.9). We find that an additional term proportional to x^'^ is 
needed in (7.10) in order to satisfy the above three constraints. 



For the data of gi{x), we will use the SMC and EMC |^ results, both being measured 



at the mean value of Qq = 10 GeV^. Following the SMC analysis we have used the new F2{x) 
structure function measured by NMC ||117|| , which has a better accuracy at low x, to update 
the EMC data. Assuming that = {Qq) = 10 GeV^ for each x bin of the gi{x, Q'^) data, a 
best least fit to gi{x) at x > 0.2 by (7.10) is found to be P3 



AM^(x,g^) = x-°-^^^(l -x)^-^^(0.928 + 0.149v^- 1.141a; + 11.612s^-^), 
A4(x,Qg) = x-°-^^^(l-x)^-^^(-0.038 - 0.43v^- 5.260x + 8.443x^-^), (7.11) 

which satisfies all aforementioned constraints. Note that we have evoluted qy{x,Q'^) from 
Q'^ = 4GeV^ to 10 GeV^ in order to compare with Aqy{x,Ql) and that the sign of Ad^ix) 
in our parametrization flips at xq = 0.496 (see Fig. 7). 

The NLO parametrization (7.11) for valence quark spin densities is obtained by assuming 
= (Qq) for each x bin of the gi{x,Q'^) data. However, a full NLO analysis should 
take into account the measured x dependence of at each x bin by considering the NLO 
evolution of parton spin distributions. At present, there already exist several such analyses 
62| , |101| , p.02| , |103|| . For example, by fitting some parametrizations for spin-dependent parton 



distributions to all available world data on gi^x^Q"^), Gehrmann and Stirling ||102|| obtained 
(set A) 

Am^(x,Q2) = 0.918A„x~°-^^^(l -x)^-^^(l -4.60v^+ 11.65x), 

Aci^(x,Q^) = -0.339AdX-°-^2°(l -a;)^-^^(l -3.48yx + 7.81x), (7.12) 

to NLO at = 4GeV^ where A„ = 1.3655 and = 3.8492 are normalization factors 
ensuring that the first moments of Auy{x) and Adjj{x) are 0.918 and —0.339, respectively. 
However, the a; ^ 1 behavior of the valence quark spin distributions (7.12) : Aujj{x)/uy{x) —>■ 
0.87 and Adv{x) / dy{x) —>■ —2.56 is not consistent with above-mentioned QCD constraint (the 
latter also seems to violate the positivity constraint). Three different NLO parametrizations 
of valence quark spin distributions are shown in Fig. 7. 

A comparison between the theoretical curve of xgi{x) fitted to the EMC and SMC data at 
X > 0.2 with the polarized valence quark distribution given by (7.11) is shown in Fig. 1 (see 
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Figure 7: NLO valence-quark spin distributions at = 10 GeV^ for three different parametriza- 
tions: (7.11) (thick sohd curve denoted by CLW |3|), the standard set of GRSV flOll (sohd curve) 
and set A of GS [|l02| (dotted curve). 



Sec. 2.2). The discrepancy between theory and experiment for g\[x) at small x is presumably 
accounted for by sea quarks and gluons. For sea-quark polarization, we know that its size is of 
order —0.10, but there is no any information on the size of gluon polarization. The sea-quark 
and gluon spin distributions cannot be separately determined by current experimental data 
102| , |101|| ; they are correlatively constrained by (7.8). In other words, vjhile the shapes of the 



spin- dependent valence quark distributions are fairly constrained by the data, the sea-quark 
and gluon spin densities are almost completely undetermined. In principle, measurements 
of scaling violation in gi{x,Q'^) via, for example, the derivative of (y'i(x, Q^) with respect 
to Q^, in next-generation experiments will allow an estimate of the gluon spin density and 
the overall size of gluon polarization. Of course, the data should be sufficiently accurate in 
order to study the gluon spin density. Meanwhile, it is even more important to probe AG'(x) 
independently in those hadron-hadron collision processes where gluons play a dominant role 
(see Sec. 8). 

As stressed in passing, the fact that gluons make no contributions to Fi in the gauge- 
invariant factorization scheme does not imply a vanishing gluon contribution to gi{x). Quite 
opposite to the naive sea-quark interpretation for (71(2;), if there is no sea polarization in 
the chiral-invariant scheme, then the size of the gluon spin component in a proton must 
numerically obey the relation AG{Q'^) = -{27r /as{Q'^))Aqf\Q'^) [cf. Eq.(4.37)] in order to 
perturbatively generate a negative sea-quark polarization AgP(Q^) via the anomaly mech- 
anism. In other words, even gluons do not contribute to Ff, the gluon spin can be as large 
as 2.5 for Aqf^ = -0.10 at = 10 GeV^ provided that Aq'^^ = 0. Recall that the gluon 
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polarization induced from quark's bremsstrahlung is positive [cf. Eq.(3.29)]. Recently, a 
full NLO analysis was performed in the chiral-invariant factorization scheme and it was 
claimed that the present data of gi are sufficient to determine the first moment of the gluon 
spin distribution, namely IS.G{Q^ = 1 GeV^) = 1.5 ± 0.8 and is roughly twice as large at 
= lOGeV^. Since the shape and size of the spin-dependent gluon distribution is k±- 



f actor ization independent, evidently there is a contradiction between the conclusions of 1^2 1 



and of |TOT], |T02|. Because AP//^(x) and AC^^^ (x) are available only in the MS scheme, it 



has been attempted to introduce a modification on NLO anomalous dimensions and hard 
coefficient functions to transfer from the GI scheme to the CI prescription ||100|| . In the CI 
scheme, ASci does not evolve with and this requires that 



Jo 



S,qq 



0. 



(7.13) 



However, this transformation cannot be unique since it is only subject to the constraints 
(7.13). Indeed, three different scheme changes have been constructed in [Q. As a conse- 
quence, the NLO evolution of polarized parton distributions in the CI scheme obtained in 



this manner |^ is ambiguous as it depends on the scheme of transformation. 



Finally, the interested reader is referred to ||118|| for a collection of polarized parton 
distributions up to 1995. 



8 Experimental Signatures of Parton Polarizations 

It is concluded in Sec. 7.2 that while the present experimental data put a useful constraint 
on the shape of the valence-quark spin distributions, the sea-quark and gluon spin densities 
are only loosely constrained. The question of what is the magnitude and even the sign of 
the gluon spin remains unanswered. In view of this, we shall survey the processes which can 
be used to probe the parton spin densities, especially for sea quarks and gluons. For the 
purposes of this section, we will use Aq and Ag to denote the spin components of quarks 
and antiquarks, respectively. 



valence-quark spin distribution Aq. 



x] 



• semi-inclusive DIS ||119| , |120| , |121|| Consider the semi-inclusive decays e + p e + 
(vr^, K^, K^,K^) +X with the longitudinally polarized lepton beam and proton target. The 
differential cross section for the production of a hadron in DIS of charged lepton has the 
form 



a 



dN^ 



dxdydz dz 



i=q,q 
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where y = u/E = {E ~ E')/E, z is the fraction of the parent parton's momentum carried 
by the final hadron h, is the number of hadrons produced with a value of z, and is 
the fragmentation function of a quark i into the hadron h. Therefore, 

^ - ^ « T.-lAg,i-,Q')Df{z,Q^ (8.2) 
Based on (8.2), two semi-inclusive asymmetries of interest are 



(8.3) 



where = lo dz Dfiz^Q^), and 



c/A^lf -'^^ /c/z + dNf^''-^~ /dz ■ 



It follows that the asymmetries in difference of deep inelastic vr^ and vr productions given 
by 

4Am^(x, Q^) - A(i^(x, 

for the proton target, and 



Af-""- (x, Q') = , r^'^-, T/'X^. (8-5) 



y,^+-7r-/ _ Am^(x,Q^) + A4(x,g^) 

for the deuteron target, are completely independent of the fragmentation function and can 
be used to extract the polarized valence quark densities. Likewise, for kaon production 

A;'-^-{x,Q') = ^^^j^, Af-^\x,Q') = ^^j^. (8.7) 

A measurement of the asymmetry in difference of K'^ and productions is thus very 
useful to test the large x behavior of Ad^{x)/dy{x), which is expected to approach unity in 
perturbative QCD, but to — | according to the Carlitz-Kaur model ||113|| . Theoretically, NLO 
corrections to semi- inclusive asymmetries were recently studied in ||121|| . Experimentally, 
asymmetries Ap^d -^p^'^ have been measured by SMC recently from which valence 
quark and non-strange quark spin distributions are extracted with the results: Au^ = 1.01 ± 
0.19 ± 0.14 and Ad^ = -0.57 ± 0.22 ± 0.11 . 

• and Z production [jl22|, |123[ In high -energy hadron-hadron collisions, the single- 
spin asymmetry Al defined by 

with da^^^'' denoting the inclusive cross section where one of the initial hadron beams is 
longitudinally polarized and has +(— ) helicity, is expected to vanish to all orders in strong 
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interactions unless some of the parton-parton scatterings involve parity-violating weak in- 
teractions. Therefore, a nonzero Al arises from the interference between strong and weak 
amplitudes and usually is small, of order 10~^ (for a recent analysis of parity-violating asym- 
metries, see [[124| , |125|| ). The only exception is the direct and Z productions in proton- 



proton collisions where a large Al of order 10% is expected to be seen at RHIC energies 
123|| . In the parton model, pp + X proceeds dominantly via ud — {ud — W~) 

^+ _ ^u{Xa,Ml,)d{Xk,Mlr)- M{Xa,Ml,)u{Xb,Ml,) 



A\ 



u{xa, M^)d{Xb, M^) + d{Xa, M^)u{xb, M^) 



^ d{Xa,Mlr)u{Xb,Mlr)+u{Xa,M^)d{Xb,Ml,) ' 

where Xa = y/r eyi\){y) , Xb = v^exp(— y), r = M^^/s. As y is near 1, we have Xa » Xb and 

^+ Au{Xa,M^) y^- AdjXg, M^) 

where Xa = e-^r. The (valence) quark spin distributions at large x and at Q"^ = thus 
can be determined at the kinematic limit ?/ ~ 1. 

antiquark and sea-quark distributions Aq{x), Aqs{x) Q 



semi-inclusive DIS ||119| , |128|| Assuming {z) > [z) > (z), we can neglect 



the strange-quark contribution in (8.2) and obtain 

. . _ ^n-TT _ 4A^(^) + ^d{x) + [Au{x) + Ad{x)]Df{z)/Dl\z) 

^ ' ^ 4«(x) + d{x) + [Au{x) + d{x)]Dt{z)/Df{z) ' ^ ' ^ 

and from A^"^ with the replacement u ^ d and u ^ d. Hence, the polarized non-strange 
antiquark distribution is determined provided that valence quark spin densities and the ratio 
Dl"^ (yz) / D'^ [z) are known. For other strategies, see [p.l9| , |120|| . Another possibility is to tag 



fast-moving K produced in semi-inclusive DIS to probe the strange-quark polarization 



128, 120 



• and Z production ||122| , |123|| In the other extreme kinematic limit y ^ — 1, we 
have Xa « Xb and the parity- violating asymmetries (8.9) become 



^+ Adjxg, M^) Au{xa,Ml^, 



^^As mentioned in the beginning of this section, we employ a different definition for quark spin densities 
here: A^s = q),{x) — q\{x) and Ag(x) = q^{x) — q^{x). A priori Aq can be different from Aqs if they are 
not produced from gluons. Based on the measurements of octet baryon magnetic moments in conjunction 



with the quark polarization deduced from DIS, it has been claimed in |126 that Aq w 0. In principle, a 



measurement of the correlations between the target polarization and the A and A polarizations in DIS will 



provide a way of discriminating between As and As (see e.g., |127|). Another nice method is to measure the 



single asymmetry in charmed meson production in the semi-inclusive DIS process as discussed below. 
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with Xa = \^/e. The parity- violating asymmetry in the kinematic region y 
information on Au{x) and Ad{x) at small x and large Q^, = Mu 



-1 provides 



Drell- Yan process [|129| , |13CI| , |131| , |132|| The double-spin asymmetry defined by 



A 



DY 
LL 



(8.13) 



da^ydQ^ + da^^/dQ^ 

measured in the Drell- Yan process pp + X, where do"^^^^-''-* designates the Drell- 

Yan cross section for the configuration where the incoming proton helicities are parallel 
(antiparallel), is sensitive to the sea spin densities. In the parton model, the asymmetry 
reads 

Aira"^ dxi dx2 ^ / Q"^ 



9Q'^s Jo Xi X2 



6 1- 



SX1X2 



^ e5Ag(xi, Q')Aq{x2, Q') + {1^2). (8.14) 



The sign of A^J is expected to be negative as Au{x) > and Au{x) < 0. A recent analysis 
of NLO effects in [|132|| indicates that the qq subprocess exhibits great perturbative stability, 
whereas the qG subprocess is important and contributes destructively. For a discussion of 
single-spin asymmetry in the Drell- Yan process, see ||122| . 



inclusive DIS with charged current [133, 134, 135 High energy lepton-proton scat- 



tering at large allows to probe spin effects in charged current interactions in the DIS 
process —>■ v{i')X. Consider the parity-violating DIS of unpolarized charged lepton on 
longitudinally polarized proton: + p ^ v + X. It is shown in [|133] that the single-spin 
asymmetry in this process is sensitive to Ady{x) and to antiquark/sea-quark spin densities: 
A'u(x), Ads{x) and As{x). With longitudinally polarized lepton and proton beams, the 
double asymmetries defined by 



+ da 



TT 



.15) 



have the expressions 134, 135 



A^~ 



A 



w+ 

LL 



Au{x) + Ac{x) - (1 - yf[Ad{x) + As{x)] 

u{x) + c(x) + (1 — yY[d{x) + s{x)] 
(1 - yY[Ad{x) + As{x)] - Au{x) - Ac{x) 



(1 — yY[d{x) + s{x)] + u{x) + c{x) 
These inclusive asymmetries are mainly sensitive to the u and d-quark fiavor. 



(8.16) 



• semi-inclusive DIS with charged current [|135|| Consider the charmed meson production 
in the semi-inclusive DIS: + p ^ z/(z/) + D + X. The main subprocesses are s + W'^ D 
and s + D. The single asymmetries are given by ||135|| 

As{x) +tan6'^A(J(a;) 



A 



W-,D. 
LL ' 



X) = 



<;'"(x) 



da 



W+,D 
Ti 



-|- da 
— da 



W-,D 
TT 

W+,D 
TT 



, W+,D , 1 W+,D 

a(7-|-j -|- acT-|-| 



s{x) + ta.n6'^d{x) 

As{x) +tan6'^A(i(x) 
s{x) + tanO^d^x) 



(8.17) 
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where 6c is a Cabibbo mixing angle. These asymmetries allow to extract the strange sea spin 
distributions As{x) and As{x) separately. It is of great interest to test if As{x) is identical 
to As(x). 



elastic vN scattering |42, 79, 136, 03 



38[| Assuming a negligible Ac, it was originally 
argued that the axial-vector form factor G^(g^) appearing in the matrix element {N\A^\N) 
for i/N elastic scattering is related to the quark polarization by Ga(0) = \{Au — Ad — As). 
Hence, a measurement of Ga(0) will determine As independently. Since the limit = is 
experimentally unattainable, the dependence of Gyi(g^) is usually assumed to have a dipole 
form. The up and up experiments in 1987 [|13|] indicated that As(0) = -0.15±0.09 |2|, 
It becomes clear now that what measured in up up scattering is the combination As — Ac 
rather than As itself [|137|| . First, contrary to the scale-dependent As, the quantity As— Ac is 
scale independent as it is anomaly free. It is possible that As at a relatively low scale is zero, 
but it evolves dramatically from the quark-model scale to the EMC scale QImc ~ GeV^ 
|141|| . Therefore, the previous interpretation for As(0) cannot be extrapolated to As((5emc) 
directly, but what we can say now is As((5^) = —(0.15 ± 0.09) + Ac{Q'^). Second, far below 
charm threshold, charmed quarks stop making contributions to DIS, but they still contribute 



to Ga via the triangle diagram. As stressed in [|1371| , the value of Ac defined in up scattering 
can only be interpretated as Ac in DIS well above the charm threshold. A new uN scattering 
experiment using LSND (Liquid Scintillator Neutrino Detector) is currently underway at Los 
Alamos (see e.g., ||136|| ). 



• semi-inclusive A production in DIS ||142| , |143| , |144|| Consider the semi-inclusive decay 
i+p i+A+X with a longitudinally polarized proton target. Since in the naive quark model 
the spin of the A is carried by the strange-quark's spin, it is expected that the negative strange 
sea polarization in a polarized proton will be transferred to the longitudinal A polarization in 
the current fragmentation region. In the simple parton model, the longitudinal polarization 
of the A is given by ||142 

dN^^/dz - dN^^/dz As(x, Q'^)AD^{z, Q^) 



Va{x,Z,Q , ^^AT/^^ + ^^Ai/^^ 

where AD^ = Df^ -Df^ = Df^ - D^l 
suggested in ||124|| a simple parametrization for 



six,Q^)Dt{z,Q^) ' 
Df^' -Dfl. Very little is known about D^{z) and AD^(z 



(8.18) 



'z): zDi 



0.5z(1.08-^) 



It is 
0.06(l-z)l 



In the absence of experimental data or a detailed theory, the construction of AD^{z) is nec- 
essarily ad hoc. Nevertheless, we expect that the polarization of the outgoing A is equal to 
that of the strange quark at z = 1. Beyond the non-relativistic quark model, the u and d 
quarks in the A are also polarized and will contribute to Va ||143|| . The longitudinal polar- 
ization of the A also can be produced in the target fragmentation region in deep-inelastic 



^^The value of Ga(0) is very sensitive to the dipole mass Ma- For example, it is shown in [140| that 
Ga{0) ^ -0.15 ± 0.07 is obtained for Ma = 1.049 ± 0.019 MeV, but the data also can be fitted with 
Ga{0) ^ and Ma = 1.086 ± 0.015 MeV. 
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vN, fl{e)N{N) scatterings and various underlying mechanisms for Va are discussed in 



144' 



gluon spin distribution AG{x) 

Phenomenological signatures of AG can be tested in various ways, as partly summarized 
in Table III. Instead of going through the details for each process, we will focus on those 
promising processes which have better signals, higher event rates and larger asymmetries. 
Since AG(Q^) increases logarithmically with Q^, it is conceivable that the effects of gluons 
will manifest in the polarized pp collider at the RHIC and in the ep collider at the HERA. 



prompt photon production [|129| , |15CI| , |151| , |152| , |153|| The double-spin asymmetry 



in the direct photon production at high pj_ in longitudinally polarized proton-proton colli- 
sions depends strongly on the polarization of gluons as the Compton subprocess Gq •yq 
dominates over qq 7G annihilation, reflecting by the fact that Ali^ grows with xp sd. fixed 
p±. This process thus provides a clean, direct, and unproblematic possibility for determining 
AG{x). Such experiments should be feasible in the near future at the RHIC. 



single-jet production ||154| , |151| , |155| , |130| , |15(j| ] In general the double-spin asymmetry 



A^l\ for a jet production in pp collisions with a transverse momentum p± is sensitive to the 
gluon spin density for x± = Ip^j ^ not too large. Since the polarized gluon distribution is 
large at small x, gluon-gluon scattering dominates the underlying parton-parton interaction 
subprocesses at small x^. As the jet momentum increases, quark-gluon scattering becomes 
more and more important due to the relatively fast decrease of the gluon spin distribution 
with increasing x. It is finally governed by quark-quark scattering at large x^. Therefore, a 
measurement of J^'^i^ in the jet momentum region where the spin asymmetry is dominated 
by qG or GG scattering will furnish important information on AG(x). 



hadronic heavy- quark production ||164| , |162| , |165| , |147|| Since the dominant subprocess 



for hadronic heavy-quark production in pp collisions is GG — > QQ, this process depends 
quadratically on AG and is hence very sensitive to the gluon spin density; it is often consid- 
ered to be the best and most realistic test on AG. 



9 Conclusions 

The new polarized DIS experiments in recent years have confirmed the validity of the EMC 
data and the controversial conclusions that the observed value of Ff, the first moment of 
gi{x), is substantially smaller than the Ellis- Jaffe conjecture and that only a small fraction 
of the proton spin comes from the quarks. However, the proton spin problem now becomes 
less severe than before. The new world average is that AS ~ 0.30 and As ~ —0.10 at 
g2 = lOGeVl The Bjorken sum rule has been tested to an accuracy of 10% level. Some 
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Table III. Various processes which are sensitive to the gluon spin distribution. 



Process 


Dominant subprocess 


References 


cn a Tin r\v J 1 1 n 1 PFtf OTiVoH 1 1 ci'M~\T\ 
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t,hrpp-ipt Drodnrtion 
PP ^ 3 jets + X 


^7^7 — > aaa 
^ qqG 


1 


155] 


fonr-ipt DroHnrtion 
— > 4 jets + X 
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1 


15C] 


two-photon production 
PP ^ 77 + X 


aa — > 
NLO: see |T6ll] 


LO: 1" 
NLO: 
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heavy quark production 
pp^ QQ + X 
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NLO: see fie^] 


LO: [" 
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charmonium production 
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1 
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main conclusions are: 

1) . There are two /c^-factorization schemes of interest: the chiral-invariant scheme in 
which the ultraviolet cutoff on the quark spin distributions respects chiral symmetry and 
gauge invariance but not the axial anomaly, and the gauge-invariant scheme in which the 
ultraviolet regulator satisfies gauge symmetry and the axial anomaly but breaks chiral sym- 
metry. The usual improved parton model calculation corresponds to the chiral-invariant 
factorization scheme. There is an anomalous gluonic contribution to due to the axial 
anomaly resided in the box diagram of photon-gluon scattering a,t kj_ — [{1 — x)/Ax\Q'^ with 
x — > 0. As a consequence, AEci is not necessarily small and Asci is not necessarily large. 
For AG ~ 2.5 at = iq QeV^ one has AEqi ~ 0.60 and Asci ~ 0. 

2) . The OPE approach corresponds to the gauge- invariant factorization scheme. Hard 
gluons do not contribute to Fi because the axial anomaly is shifted from the hard photon- 
gluon cross section to the spin-dependent quark distribution. However, it by no means 
implies that AG vanishes in a polarized proton. AEgi is small because of the negative 
helicities of sea quarks. The chiral-invariant and gauge-invariant factorization schemes are 
explicitly shown to be equivalent up to NLO since Ag's and Acr^^^'s in these two different 
schemes are related by (4.36) and (4.33), respectively. As far as the first moment of gi{x) is 
concerned, the anomalous gluon and sea-quark interpretations are thus on the same footing. 

3) . Contrary to the gauge-invariant Aq'gi, AG and chiral-invariant Agci cannot be ex- 
pressed as matrix elements of local and gauge-invariant operators. Nevertheless, gauge vari- 
ant local operator definitions do exist; in the light-front coordinate and in the light-front 
gauge = (or in the infinite momentum frame and in temporal axial gauge), AG has 
a local operator definition given by (4.55) or (4.56), and Agci by (4.57). By contrast, they 
can be also recast as matrix elements of string-like gauge-invariant but non-local operators. 

4) . The U(l) Goldberger-Treiman relation (5.6) in terms of the r]Q remains totally un- 
changed no matter how one varies the quark masses and the axial anomaly, while its two- 
component expression (5.17) is identified with the connected and disconnected insertions [see 
(5.21)]. We have determined the physical coupling constants g^,^j^ and g^j^^^ from the GT 
relations for g'^ and g\ and found that g^,^j^ = 3.4 and g^j^^^ = 4.7 . 

5) . For massless sea quarks there are two mechanisms allowing for quark helicity fiip and 
producing sea-quark polarization: the nonperturbative mechanism due to instanton-induced 
interactions (see Sec. 4.4) and the perturbative way via the axial anomaly (see Sec. 4.3). The 
sign of the sea-quark helicity generated by hard gluons via the latter mechanism is predictable 
in the framework of perturbative QCD: It is negative if the gluon spin component AG is 
positive. The lattice calculation indicates that sea polarization is almost independent of 
light quark fiavors; this empirical SU(3)-fiavor symmetry implies that it is indeed the axial 
anomaly, which is independent of light quark masses, that accounts for the bulk of the 
helicity contribution of sea quarks. 

6) . A full and consistent next-to- leading order analysis of the gi{x, Q^) data just became 
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possible recently. We have to work entirely in the gauge- invariant factorization scheme for the 
NLO analysis since the NLO polarized splitting functions are available only in this scheme. 
While the shapes of the spin-dependent valence quark distributions are fairly constrained by 
the data, the sea-quark and gluon spin densities are almost completely undetermined. It is 
thus very important to probe AG{x) independently in the hadron-hadron collision processes 
where gluons play a dominant role. The most promising processes are: prompt photon 
production, single-jet production and hadronic heavy-quark production in pp collisions. 

7). As for the spin sum rule | = + + LI + L'^, the only spin content which 
is for sure at present is the observed value AE ~ 0.30 at Q"^ — lOGeV^. The relativistic 
quark model predicts that AE = 0.75 and L| = 0.125. Recent lattice calculations imply 
that relativistic quark model results are recovered in the valence approximation. The quark- 
model's value of 0.75 for AE is reduced to the "canonical" value of ~ 0.60 by negative 
spin components of cloud quarks, and reduced further to ~ 0.30 by the negative sea-quark 
polarization. That is, the deviation of from unity expected from the non-relativistic quark 
model is ascribed to the negative spin components of cloud and sea quarks and to relativistic 
effects. The "valence" contribution as conventionally referred to is actually a combination 
of cloud-quark and truly valence-quark components. It is thus important to estimate the 
cloud-quark polarization to see if it is negative in sign and comparable in magnitude to the 
(one-fiavor) sea-quark helicity. In the asymptotic limit, J^ioo) = |AE(cxo) -|-L^(oo) ~ ^ and 
J^(oo) = AG(oo) -|- L^((X)) ~ ^. If the evolution of J| and is very slow, we will have 
L^(10 GeV^) ~ 0.10, which is close to the quark-model expectation. The growth of AG with 

is compensated by the gluon orbital angular momentum, which also increases like In 
but with opposite sign. 
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